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ABSTRACT 

We present the results from a survey of 57 low-redshift Abell galaxy clusters to study the radial 
dependence of the luminosity function (LF). The dynamical radius of each cluster, r2oo> was estimated 
from the photometric measurement of cluster richness, B gc . The shape of the LFs are found to correlate 
with radius such that the faint-end slope, a, is generally steeper on the cluster outskirts. The sum 
of two Schechter functions provides a more adequate fit to the composite LFs than a single Schechter 
function. LFs based on the selection of red and blue galaxies are bimodal in appearance. The red LFs 
are generally flat for —22 < Mr c < —18, with a radius-dependent steepening of a for Mr c > —18. 
The blue LFs contain a larger contribution from faint galaxies than the red LFs. The blue LFs have a 
rising faint-end component (a ~ —1.7) for M# c > —21, with a weaker dependence on radius than the 
red LFs. The dispersion of M* was determined to be 0.31 mag, which is comparable to the median 
measurement uncertainty of 0.38 mag. This suggests that the bright-end of the LF is universal in 
shape at the 0.3 mag level. We find that M* is not correlated with cluster richness when using a 
common dynamical radius. Also, we find that M* is weakly correlated with BM-type such that later 
BM-type clusters have a brighter M*. A correlation between M* and radius was found for the red 
and blue galaxies such that M* fades towards the cluster center. 



Subject headings: Galaxies: clusters: general 
— Galaxies: evolution 



Galaxies: luminosity function — Galaxies: formation 



1. INTRODUCTION 

The study of the formation and evolution of galaxies 
is a fundamental avenue of research in the process of un- 
derstanding astrophysical and cosmological issues. How 
galaxies form and evolve can be studied using a variety 
of techniques, one of those being the galaxy luminos- 
ity function (LF). The galaxy LF, assuming that galaxy 
mass-to-light ratios are nearly constant for similar types 
of galaxies, can potentially provide a direct link to the 
initial mass function and hence the distribution of den- 
sity pertanbationsjhatjre thought to give rise to galax- 
ies (jPress fe Schecterlll974[ ) . Since most galaxies are not 
isolated entities, evolutionary processes, in addition to 
those expected for an aging stellar population, can occur 
as galaxies interact with their environment. 

The galaxy LF — the number of galaxies per unit vol- 
ume in the luminosity interval L to L+dL — can be used 
as a diagnostic tool to search for changes in the galaxy 
population. In particular, the LF for cluster galaxies can 
help ascertain the influence of the cluster environment 
on the galaxy population. For example, a change in the 
shape of the LF with respect to cluster-centric radius 
provides important insight into the dynamical processes 
at work in the cluster environment. 

A central theme in the early studies of the galaxy clus- 
ter LF has bee n to determine whether the LF is universa l 
in shape (e.g., iHubbld HIM libel [l962l lOemlerl fl97l . 
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While introducing the mode rn form o f the L F, the so- 
called "Schechter Function" , ISchechterl (|1976f ) suggested 
that the cluster LF is universal in shape, and can be 
characterized with a turnover of Mg — —20.6 + 5 log hso 
and a faint-end slope of a = —1.25. Further support for a 
universal LF has been provid e d by several studies such as 
iLuggerl (fl98 l. iCollessI (Il989h lOaidosI (fl997ri . lYaeTeTaTl 
(2002), and iDe Propris et alj (2003?) . In contrast, sev- 
eral studies have show n that the shape of the cluster 
LF is not universal (e.g.. [Godwin fe Peachlll977HDresslerl 
19781: lLopez-Cruz et al.1ll997t iPiranomonte et all 120011: 
Hansen et al.l 120051 : iPopesso et al.l |2006|) . One expects 



that the LF depends on cluster-centric radius since the 
mixture of galaxy morphological types should vary with 
radius, as imp lied by the morphology-density relation 
(|Dresslerlll980l) . Since different morphological type s are 
characterized by different LFs (Bing geii et al.lfl988l ) the 
cluster LF (integrated over all galaxy types) should not 
be universal. Indeed, some studies have provided evi- 
dence that the cluster LF does vary w i th cluster-centric 
radius (e.g..lBeiiersbergen et al.l' 2002a: G oto et aT1l2005l : 
lHansen et al.ll2005l : IPopesso et al.ll2006 ). 

The main goal of this paper is to investigate the change 
in the cluster i? c -band LF as a function of cluster-centric 
radius. To avoid the inherent bias that has plagued nu- 
merous studies, the cluster LF will be compared based 
on scaling relative to the dynamical radius, r2oo- The 
use of a dynamical radius for comparing galaxy pop- 
ulations for a sample of clusters, provides one of the 
most robust, least-biased, photometric survey yet pub- 
lished on the LF of Abell clusters. Directly compar- 
ing LFs that sample only the cluster core region with 
other cluster LFs that extend to the outskirts, will suf- 
fer from radial sampling bias given that the shape of the 
LF h as been shown to depend on c luster-centric radius 
(e.g., iChristlein fe ZabludofB 120031 lHansen et all [20051: 
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iPopesso et all [2006). A direct comparison of the galaxy 
population with respect to cluster-centric radius based 
on ^200, will help to accurately measure the change in 
the properties of cluster galaxies as a function of global 
environment. These data will also provide information to 
help settle the long-standing debate regarding the univer- 
sality of the cluster galaxy LF and the properties of the 
faint dwarf galaxy component. 

This paper is the third in a series resulting from a 
multi-color imaging survey of low-redshift (0.02 < z < 
0.2) Abell clusters. The paper is organized as follows. In 
§2 we present a brief overview of the sample selection and 
data reduction procedure. In §3 we describe the method- 
ology for generating the galaxy cluster LF. In §4 we ex- 
amine the LF for different color-selected galaxy popula- 
tions of our cluster sample. Discussion and conclusions 
are presented in §5. Finally various selection effects 
and biases are explored in the Appendix. Further de- 
tails regarding sample selection, observations, image pre- 
processing, cata^gSj_and^iiding_charts can be found in 
lLopez-Cruzl ((1997), Ba rkhouse] (|2003f) . and Barkhouse et 
al. (2007a; in preparation, Paper I of this series). A de- 
tailed discussion of the color-magnitude relation (CMR) 
of early- t ype galaxies using this survey can be found in 
Paper II (|L6pez-Cruz et aL| [2004) . Paper IV characteris- 
tics the cluster galaxy luminosity and color distribution 
by examining the dwarf-to-giant ratio and the galaxy red- 
to-blue count ratio (Barkhouse et al. 2007b, in prepara- 
tion). Recent observations suggest that the best cosmo- 
logical model is characterized by — 0.3, £l\ ~ 0.7, 
and H ~ 70 km s" 1 Mpc -1 (e.g., iSpergel et al.|[2003h . 
Since the effects of curvature and dark energy are negligi- 
ble at low-redshifts (z < 0.2) and to allow direct compar- 
isons with previous studies, we have set for convenience, 
unless otherwise indicated, Hq = 50 km s -1 Mpc -1 and 
qo = throughout this paper. 

2. OBSERVATIONS AND DATA REDUCTIONS 

The galaxy cluster sample utilized for this paper is 
identical to that descri bed in Paper II of this series 
(|L6pez-Cruz et all 12004). We summarize the observa- 
tions and data reductions below. 

The galaxy cluster sample is composed of Abell clus- 
ters selected ma i nly f rom the X-ray compilation of 
Uones fc Formanl (|1999t) . The primary cluster sample 
was selected based on the following criteria; (1) clusters 
should be at high galactic latitude, \b\ > 30°; (2) their 
rcdshifts should lie within the range 0.04 < z < 0.20; 
(3) the Abell richness class (ARC) should, preferably, be 
> 0; and (4) the declination S > -20°. Some ARC = 
clusters were included in the final sample due to the 
lack of suitable clusters at certain right ascensions dur- 
ing the observations. This sample includes 47 clusters of 
galaxies observed in B 7 and Kron-Cousins R c and I at 
KPNO with the 0.9 m telescope usi ng the 2048 x 2048 
pixel T2KA CCD (0.68" pixel" 1 ) ([Lopez-Cruzl Il997t 
lYee fc L6pez-CruzHl999l lLopez-Cruzl 1200 ll hereafter the 
LOCOS sample; Lopez-Cruz et al. 200 7, in preparation) . 

A sub-sample of eight clusters from iBarkhousel ([2003) 
is included to complement our 47-cluster sample by cov- 
ering the low-redshift interval from 0.02 < z < 0.04. 
These data were obtained at KPNO with the 0.9 m tele- 
scope using the 8K MOSAIC camera (8192 x 8192 pix- 



els; 0.423" pixel - ). The clusters for this sample were 
selected using the previous criteria except that ARC = 
clusters were not preferentially excluded. In a ddition, 
two cl usters imaged in B and R c are included from lBrownl 
using the same instrumental setup as the LO- 
COS sample and selection criteria as the MOSAIC data. 
All clusters in the sample are detect ed in X-rays and 
are fo und to have a prominent CMR (|L6pez-Cruz et al] 
l20f)l . 

The integration times for our 57-cluster sample varies 
from 250 to 9900 s, depending on the filter and the red- 
shift of the cluster (only the B- and i? c -band data are 
considered here). Control fields are also an integral part 
of this survey. For this study we use a total of six con- 
trol fields in both the B and R c filters. The control fields 
were chosen at random positions on the sky at least 5° 
away from the clusters in our sample. These control fields 
were observed using the MOSAIC camera to a compara- 
ble depth and reduced in the same manner as the cluster 
data. All observations included in this study were carried 
out during 1992-1993 and 1996-1998. 

Processing of the 8k mosaic images were done using 
the mscred package within the IRAF environment. The 
photometric reduction was carried out using t he program 
PPP (Picture Processing Package; lYedll99ll ). which in- 
cludes algorithms for performing automatic object find- 
ing, star/galaxy classification, and total magnitude de- 
termination. A ser ies of improvements to PPP described 
in lYee et all (|1996l ) that decreases the detection of false 
objects and allows star/galaxy classification in images 
with a variable point-spread-function (PSF) was utilized. 

The object list for each cluster is compiled from the 
R c frames. The R c frames are chosen because they are 
deeper than the images from the other filters. Galaxy 
total magnitudes are measured with PPP using a curve- 
of-growth analysis. The maximum aperture size ranged 
from 20" for faint galaxies (R c > 18.5) to as large as 
120" for cD galaxies in z ~ 0.02 clusters. An optimal 
aperture size for each object is determined based on the 
shape of t he curve-of-growth using criteria described in 
lYed (fl991). The photometry of galaxies near the cluster 
core was carried out after the cD and bright early-type 
galaxies had been r emoved using profile-modeling tech- 
niques developed bv iBrownl (|1997l ). 

Galaxy colors were determined using fixed apertures 
of 11.0/i^Q kpc on the images of each filter at the red- 
shift of the cluster, sampling identical regions of galaxies 
in different filters, while imposing a minimum color aper- 
ture of ~ 3 times the full-width-half-maximum (FWHM) 
in order to avoid seeing effects (average seeing ~ 1.5"). 
The overall internal accuracy in the color determinations 
is ~ 0.005 mag in B — R c for bright objects. The errors 
for faint objects can be as large as 0.5 mag in B — R c . We 
note that the total magnitude of a galaxy is determined 
using the growth curve from the R c image, while the to- 
tal magnitude in the B image is determined using the 
color differ ence with respect to the R c image (for more 
details, see lYe^fl99lh . 

Star/galaxy classification was performed within PPP 
using a classifier that is based on the comparison of the 
growth curve of a given object to that of a reference PSF. 
The reference PSF is generated as the average of the 
growth curves of high signal-to- noise ratio (S/N), non- 
saturated stars within the frame. The classifier measures 
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the "compactness" of an object by effectively comparing 
the ratio of the fluxes of inner and outer parts of an 
object with respect to the reference PSF. 

Instrumental magnitudes are calibrated to t he Kron- 
Cousins system by observing standard stars from lLandoltl 
(1992). Due to the large field, up to 45 standard stars 
can be accommodated in a single frame. The color prop- 
erties of the standard stars cover a large color range that 
encompasses those of elliptical and spiral galaxies. The 
standard stars are measured using a fixed aperture of 30 
pixels for the LOCOS frames and 32 pixels for the MO- 
SAIC data. These aperture sizes are selected as being 
the most stable after measuring the magnitudes using a 
series of diameters. We adopt the average extinction co- 
efficients for KPNO and fit for the zero points and color 
terms. The rms in the residuals of individual fittings is 
in the range 0.020 — 0.040 mag, which is comparable to 
the night-to-night scatter in the zero points. This can 
be considered as the systematic calibration uncertainty 
of the data. 

The final galaxy catalogs were generated using the in- 
formation and corrections derived previously. For data 
obtained under non-photometric conditions, single clus- 
ter images were obtained during photometric nights in 
order to calibrate the photometry (three clusters in to- 
tal). The completeness limit for each field is based on 
a fiducial 5er limit determined by calculating the magni- 
tude of a stellar object with a brightness equivalent to 
having a S/N = 5 in an aperture of 2". This is done by 
scaling a bright unsaturated star in the field to the 5ct 
level. Since the 5<r limit is fainter than the peak of the 
galaxy count curve, and hence below the 100% complete- 
ness limit for galaxies, a conservative 100% completeness 
limit is in general reac hed at 0.6 — 1.0 mag brighter than 
the 5cr detection. See lYed f|1991h for a detailed discus- 
sion of the completeness limit relative to the 5a detection 
limit. 

Galaxy colors and magnitudes were corrected for the 
extinction produced by our Galaxy. The values of 
the Galactic extinction coeffi cients were calculated from 
the iBurstein fc Heilesl (|1982h maps using the reported 
E(B — V) values, or directly from the Ab tabulations 
for bright galaxies (jBurstein fc Heilesl Il984h with coor- 
dinates in the vicinity of our pointed observations using 
NED. Extinction values used for each cluster are pro- 
vided in Paper I. 

3. THE LUMINOSITY FUNCTION 

3.1. Fitting the Schechter Function 

The galaxy cluster LF in modern times has mainly 
been paramete rized using the Schechter function 
([Schechterlll976h . This function has the form 

<fi(L) dL = (/>* {L/L*) a exp(-L/L*) d(L/L*), (1) 

where 4>(L) dL is the number of galaxies per unit vol- 
ume in the luminosity interval L to L + dL, <p* is the 
number per unit volume, and L* is the "characteristic" 
luminosity. This function is characterized by having an 
exponential shape at the bright-end and a power-law like 
feature, whose slope is measured by a, at the faint-end 
(for example, see Figure 1). By introducing the change 
of variables M - M* = -2.5 log(L/L*), equation 1 can 



be written in terms of absolute magnitude as 

n(M) dM = kN * e lH*+i)(M*~M)~ CX p{k(Ar-M)}] dM ^ 

(2) 

where k = 0.4 In 10 (cf.. lCollesdfl989h . 

The Schechter function is fit to the clus t er ga laxy 
counts following the procedure in lLopez-Cruzl (|1997h . In 
summary, the function parameters M* , N*, and a are 
estimated by performing a x 2 minimization of the form 

x2 = E w^3! (3) 

i 1 

where iV,; is the net galaxy counts in the i th bin of the 
observed LF, Nf is the expected number of counts in 
the i th bin of width AM, and o~i is the uncertainty of 
the counts in the i th data bin. The expected number of 
counts in the i th bin, corrected for the loss of informat ion 
due to the finite bin size, is given by (|Schechter!ll9"76l) 

Nf = n(Mi) AM + n- (M,) AM 3 /24, (4) 

where the derivative is with respect to absolute magni- 
tude, and AM is the bin width. This correction is de- 
rived by Taylor-expanding n(M) to a third order about 

the bin's center and integrating Nf = n(M)dM . 

The uncertainty, a*, is taken to be (cf.. lLuggerlll986f ) 

n = [(Nf + Nbi) + 1.69JV&] 1/2 , (5) 

where Nu is the background counts in the i th bin, 
{Nf + Nbi) 1 / 2 is the Poisson uncertainty in the uncor- 
rected LF counts, and the second term in the square 
root expression is the measured field-to-field variation 
per bin in the background field counts (see §3.2). The 
uncertainty of the observed net galaxy counts due to cos- 
mic variance is thus taken into account when fitting the 
Schechter function. For a single Schechter function fit, 
TV* is fixed so that the total number of observed galax- 
ies in the data set is equal to the number predicted by 
the Schechter function. The parameters M* and a are 
obtained by minimiz ing equation 3 us ing the Levenberg- 
Marquardt method ([Press et alJll992f ). 

3.2. Background Galaxy Correction 

The i? c -band LF is constructed statistically by sub- 
tracting a background galaxy population from the clus- 
ter galaxy counts. This method, in contrast to measuring 
the redshift of individual galaxies in the cluster field, re- 
lies on an accurate determination of the background field 
population. This statistical ap proach has been used in 
numerous studies to date (e.g., [ Ocmlcr 1974: Sch echterl 
[1975 ICollessI [19891: JPrTver et all [1998L lYagi et alJ [20021 
lAndreon et al.ll2004f) . The modal field-to-field variation 
per magnitude bin has been measured to be ~ 30% above 
Poisson statistics among the six background fields. We 
use this value to approximately account for the additional 
uncertainty in the background counts due to field-to-field 
variations in equation 5. Although the use of galaxy red- 
shifts would provide a more robust determination of the 
cluster LF, the relatively modest variation of the back- 
ground counts, with respect to the cluster counts, makes 
the statistically derived cluster LF valid. A similar con- 
clusion has al so been reached by a number of independent 
studies (e.g.. lDriver et al.lll998l ). 
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Several studies have examined the effect on the de- 
rived cluster LF using a global background galaxy field 
corre ction versus one measured locally for each cluster 
(e.g.,lGoto et al.ll2002l: lHansen et aLll2005t IPopesso et al. 
120051: iGonzalez et al.ll2006f ). For example. iPopesso et al. 
(|2005h has shown for a study of 69 clusters based on 
SDSS DR2 data, that there is no significant difference 
in the measured cluster LF using either a global or local 
background subtraction technique. 

Contamination from the 2-d projection on the sky of 
a distant cluster in the field-of-view of the target cluster 
can prove to be problematic by skewing the LF, especially 
at the faint-end where galaxies from the background 
cluster directly add to the desired LF. Fortunately, the 
CMR can help minimize this contamination by identi- 
fying the early-type red sequence of the target cluster 
(|L6pez-Cruz et alll2004D . The effect of background clus- 
ters on the desired LF can thus be reduced by selecting 
an appropriate color cut, thus eliminating objects red- 
der than the cluster red sequence. This method can also 
help to locate foreground clusters which, given the low 
rcdshift of our cluster sample, are not a significant con- 
cern for this study. 

To minimize contamination from background galaxies, 
we cull galaxies that are 0.22 mag redward of the CMR 
(i.e., > 3.0 times the average B — R c dispersion of the 
cluster red sequence). The dispersion of the cluster red 
sequen ces are tabulated in Table 1 of lLopez-Cruz et all 
(2004) and histogram representations of the rectified 
B — R c color distri butions are presented in Paper I (see 
also Figure 1 from lLopez-Cruz et al.ll2004D . In several 
cases (e.g., A2152), a color cut < 3er is used if a second 
red sequence from a more distant cluster (and hence at 
a redder B ~ R c ) is apparent in the color-magnitude dia- 
gram (CMD) of the target cluster. For faint magnitudes, 
we cull galaxies redder than 2.5o-b~r c if 2.5 times the 
average uncertainty in the galaxy B — R c is redder than 
3.0 times the dispersion of the cluster red sequence. 

3.3. K- Correction 

K-corrections are applied using a single parametriza- 
tion based on early- type galaxies, which dominant the 
cluster galaxy population Dressier (198QJ). At the low 
rcdshift of our sample, the difference between early- and 
late-type galaxies is minimal for the i? c -band. In gen- 
eral, the maximum K-correction applied was ~ 0.2 mag. 
The K-correction adopted for each cluster is tabulated 
in Paper I. All LFs presented in this paper have been 
extinction- and K-corrected, and no attempt has been 
made to correct for individual internal galaxy absorp- 
tion. 

3.4. Cluster Dynamical Radius 

Studies of the properties of galaxy clusters, such as 
the CMR and LF, have been routinely compared on a 
cluster-by-cluster basis. Nearly all of these studies de- 
fine a "cluster" based on the total area covered by the 
telescope detector (e.g. , iDresslerll 1 9 78t iLopez-Crudl 1 99 7l) 
or by using a sp ecific physical length (e.g., 1 h±Q Mpc; 
lYagi et aTll2002f). The lar ge variation in cluster richness 
( Yee fe Lopez- Cruzl 1 19991 ) inhibits the usefulness of the 
above techniques for robustly comparing cluster prop- 
erties. Some authors have attempted to "normalize" 



clusters by directly comparing only clusters of compa- 
rable class (e.g., Bautz-Morgan typ e) or by weighting 
each cluster according to richness dGarilli et alj 119991 : 
iPiranomonte et alj 120011: iDe Propris et al.ll2003fh The 
use of a variety of different methods have certainly con- 
tributed to conflicting results that have emerged from 
past investigations reg arding measurements such as the 
universality of the LF (|Luggerill986l ; iDriver et alj|1998t 
IPopesso et alll2006t l. 

As a means of computing a "dynamical" radius within 
which cluster characteristics can be robustly compared, 
the r2oo radius was calculated for each cluster. The r 2 oo 
radius marks the size of a cluster within which the av- 
erage density is 200 times the cr itical density, and fol - 
low s from the defini t ion u sed in ICarlberg et al.l (|1997f ) 
and I Yee fc Ellingsonl (|2003h . The r 20 o radius is expected 
to co ntain the bulk of th e virialized mass of a cluster 
(e.g., ICole fe L ace vl [19961) and is used in this study as 
a scaling factor to compare cluster features. The use 
of this type of "normalization" allows us to compare 
galaxy cluster populations in a less biased fashion, es- 
pecially those properties which are a function of cluster 
richness and cluster-centric ra dius. This approach has 
recently been i mplem ented by Hans en et alj (|2005f ) and 
IPopesso et al.1 (|2006|) to study the properties of a sample 
of clusters/groups from the SDSS. 

The procedure for determining r 2 oo involves the calcu- 
lation of the velocity dispersion, which requires redshift 
measurements for a number of cluster galaxies. The data 
available for this study do not include redshift informa- 
tion for cluster members, thus an alternative method was 
employed to estimate T2oo for each cluster. 5 This method 
relies on the correlation betwee n B ac and r 2 oo as mea - 
sured for the CNOC1 sample (|Yee fe Ellingsonl l2003h . 
The B gc parameter is a measure of the cluster-center 
galaxy correlation amplitude, and ha s been shown to be 
a robust estimator of cluster richness (|Yee fc Lopez-Cruzl 
119991: lYee fc Ellingsonl |2003[ ) . The measured values of 
B gc for our cluste r sample are calculated using the 
method outlined in I Yee fc Lopez-Cruzl {1999) . Figure [2] 
depicts the relationship betwee n rim and B nr fo r 15 clus- 
ters from the CNOC1 surve y (|Yee et al.lll996D . adopted 
from I Yee fc Ellingsonl (|2003h . We note that the r 20 o vs. 
Bgc figure in Yee & Ellingson (se e their Figure 5) used 
the less well-determined r' 200 from ICarlberg et al.l (|1997l 
see their Table 1). This explains the decrease in the 
scatter of r 2 oo vs. B gc for our Figure [2] compared to the 
corresponding figure in Yee & Ellingson. A fit to these 
data yields 

log r 200 = (0.48 ± 0.10) log B gc - (1.10 ± 0.31), (6) 

where r 2 oo has units of Mpc, and B gc units of Mpc 1 ' 8 . 
The rms scatter in the derived values of r 2 oo is on the or- 
der of 15%. The fit was performed using the bisector bi- 
variate c orrelated errors and intrins ic scatter (BCES) es- 
timator (jAkritas fc Bershadvl 119961 ) since this algorithm 
accounts for uncertainties in both variables. The clus- 
ter MS 1455+22 (open circle in Figure [2]) was excluded 
from the fit since it is ~ 3cr from the expected relation 
(see Yee & Ellingson 2003 for a detailed discussion of MS 
1455+22). 

Only approximately 30% of the 57 clusters have published ro- 
bust velocity dispersions. 
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The estimation of r 2 oo for each cluster is accomplished 
by calculating the value of B gc directly from the galaxy 
cluster photometric catalog and then applying equation 
6 to determine r 2 oo- Table 1 lists the values of B gc and 
^200 for each cluster used in this study. The tabulated 
uncertainty in the value of r 2 QQ is calculated directly from 
t he 15% rms scatter 

lYee fc Lopez-Cru^ (|1999D showed that a subset of our 
Abell cluster sample has a similar relationship between 
cluster velocity dispersion (a v ) and B gc as the CNOC1 
sample, and since r2oo is estimated from a v , the r 2 $Q-B gc 
relation for the Abell clusters should be similar to that 
of the CNOC1 sample. 

To test the validity of our results, we compare the val- 
ues of r2oo fr o m Ta bl e 1 for those c l usters i n common with 
iRines et al.l (I2003D. iMiller et all (|2OO50 . iPopesso et all 
(|2007ft . and lAguerri et alj (|2007t ). For the 3 clusters in 
common with Rines et al., we find a mean difference in 
r 2 oo of —0.17 ±0.17 Mpc (all physical length scales have 
been converted to our distance scale), where our values 
are greater on average (the rms is given as the uncer- 
tainty). For the 10 clusters in common with Miller et al., 
we find a mean difference of 0.73 ± 0.73 Mpc. The mean 
difference for the 9 clusters in common with Popesso et 
al. is 0.25±0.32 Mpc, while the average difference for the 
12 clusters in common with Aguerri et al. is —0.26 ±0.62 
Mpc. The larger discrepancy with the Miller et al. sam- 
ple is due to the uncommonly large r2oo values (r 2 oo > 4 
Mpc) for the three richest clusters. In fact, Miller et 
al. cautions that the radius within which the density 
measurements have been made to determine r 2 Qo may 
be inaccurate. We believe that the Miller et al. r 2 oo 
values are biased-high for the more massive clusters in 
our comparison sample. If we restrict our analysis to the 
combined Rines et al., Popesso et al. and Aguerri et al. 
samples, we find a mean r2oo difference of —0.07 ± 0.53 
Mpc. Thus, the values of t 2 qq derived via equation 7 
are reasonable. (See the discussion in §A.5 regarding the 
effect on the derived LF for a 15% scatter in r 2 Qo-) 

4. RESULTS 

4.1. Individual Cluster Luminosity Functions 

The i? c -band LFs for the 57 clusters presented in this 
paper are depicted in Figure 1 for galaxies brighter than 
the 100% completeness limit. To help facilitate the com- 
parison between clusters, the LFs are generated using 
galaxy counts within a radius of (r/r2oo) < 0.4 from 
each cluster center. The cluster center is normally se- 
lected using the brightest cluster galaxy (BCG) or, when 
some doubt exists, the brightest early-type galaxy that 
is closest to the X-ray cen troid as given, for example, 
by iJones &: Formanl (fl999) . Figure 1 includes the R c - 
band LF for A496 and A1142 from data obtained by 
iBrownl (|1997t ). The small size of the detector (2kx2k) 
and the low redshift, limits the LF coverage to within a 
cluster-centric radius of (r/r2oo) < 0.2 for Abell 496 and 
(r/r 200 ) < 0.3 for Abell 1142. 

The cluster LFs presented in Figure 1 have each been 
fit with a Schechter function in the range —24 < Mr c < 
—20, with the faint-end slope fixed at a = —1. This 
has been done to help serve as a reference point for 
comparing individual clusters. The fitted value of M* 
for each cluster is tabulated in Table 1 and the best- 
fit Schechter function is represented by the solid line in 



Figure 1. The fitting of the LFs does not include the 
BCG. The presence of these galaxies is easily noticed 
by their affect on the brightest magnitude bin, whose 
value is usually offset from the best-fit Schechter func- 
tion. ISchechterl (|1976h remarked that BCGs do not seem 
to be a natural extension of the cluster LF (see also, 
iSandaed [1971 iDresslerl fl97l lLoh fe Strauss! l200l . A 
fit to his composite LF for a sample of 13 clusters was 
more robust when the BCG was excluded from each clus- 
ter. This has led to numerous debates on the forma- 
tion mechanism of BCGs (e.g.. iGeller &: Peebles! 119761 : 
lBhavsarlfl989l iBernstein fc Bhavsarll2001h . Since most 
studies of cluster LFs exclude the BCG from the LF fit, 
BCGs will not be included in subsequent LF analysis un- 
less otherwise noted (see additional discussion in §A.6). 

Visual inspection of the individual cluster LFs in Fig- 
ure 1 indicates that, in general, the bright-end appears 
to be well-fit by a Schechter function, with a rising faint- 
end (Mr o > —19) of various strength. The presence of 
a "flat" LF (a — —1.0) at the faint-end usually occurs 
for clusters in which the uncertainty in the net galaxy 
counts at the faint-end is large, or the LF has not been 
sampled to a sufficient depth to reveal a rising faint-end 
component, although we can not rule out intrinsically 
flat LF clusters. 

An additional feature, visible for several cluster LFs, is 
a "dip" in the galaxy counts at -20 < M Rc < -19. This 
characteristic is most-prominent for A84, A154, A634, 
A690, A1291, A1569, A1656, A1795, A2384, and A2556. 
A possible cause of this feature is the variation in the 
ratio of galaxy types that comprise the individual cluster 
galaxy population. This may result from the fact that 
the elliptical and spiral galaxy LF is better described 
using a Gaussian function while the dwarf gal axy LF is 
"Schechter-like" in shape (Binggcli et al. 1988|). 

4.2. Composite Luminosity Function 
4.2.1. Total Luminosity Function 

The statistical subtraction of a background galaxy 
population can drastically affect the accuracy of deter- 
mining the shape of the cluster LF when the net number 
of cluster galaxies is small. Generally, this will be an 
important factor for galaxies located in the outskirts of 
clusters. To reduce the uncertainty in the shape of the 
LF at large cluster-centric radii, cluster galaxy counts 
have been combined to form a composite LF. This also 
averages out any apparent variations in the shape of the 
individual LFs due to cosmic variance in the background 
counts. To provide an adequate coverage of the faint- 
end of the LF, we have selected a sub-sample of 29 clus- 
ters for the composite LF that are 100% complete to 
M Rc = —16.5 (absolute i? c -band completeness limits, 
Mft° m , are tabulated in Table 1). Following the proce- 
dure for combining c luster counts in adjacent magnitude 
bins (|Schechterlll976f ). clusters complete to Mr c = —17.0 
(13 additional ones) are also included in our composite 
LF. This has been accomplished by scaling the num- 
ber of net galaxy counts in the faintest magnitude bin 
(-17.0 < M Rc < -16.5) by the ratio (N 2 /Ni); where 
Ni is the total net galaxy count to M Rc = —17.0 for the 
29 clusters complete to at least M Rc = —16.5, and N 2 
is the total net galaxy count to M Ra — —17.0 for the 42 
clusters complete to at least M Rc = —17.0. By following 
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this method, we are able to construct a composite LF 
which is complete to Mr c = —16.5 and contains galaxy 
counts from 42 individual clusters. 

To search for differences in the LF which may corre- 
late with cluster properties, we present the composite 
i? c -band LF in Figure[3]for four different radial bins, cen- 
tered on the BCG, extending out to (V/r2oo) = 1. Exam- 
ination of Figure [3] clearly shows that a single Schcchtcr 
function is inadequate to fully describe the shape of the 
composite LF at any radius. The sum of two Schechter 
functions has therefore been used to model the shape of 
the LF, with the resultant fit given by the solid line in 
Figure [H 

Several recent studies have determined that the 
sum of two Schechter functions provides a more ad- 
equate fit to t he cluster LF than a single Schechter 
function (e.g., iDriyer et alj 11994 iHilker et all 120031 : 
IGonzalez etahl 120061: iPopesso et all |200j). Alternative 
LF fitting functions include a Gaussian for the bright- 
end a nd a single Schechter funct i on for the faint-end 
("e.g.. iThompson fc Gregory! Il993t iBiviano et all 119951 : 
iParolin et al.ll2003h. a single po wer-law fit to the faint- 
end (e.g.. iTrentham et all 20011). and a n Erlang plus a 
Schechter function (Biviano et al.|[l995T ). 

For Figure [3l a single Schechter function was first fit 
to the bright-end from —24 < Mr c < —20, and a second 
Schechter function from —24 < Mr c < —16.5. Due to 
the degeneracy involved in fitting two Schcchtcr func- 
tions, the faint-end slope of the first Schechter func- 
tion (for the bright-end) was fixed at ot\ = —1.0. This 
procedure is justified since a fit to all 57 clusters for 
-24 < M Rc < -20 yields a = -0.96±0.04. Limiting the 
analysis to the 42 clusters used to form the composite LF 
gives a = —0.95 ± 0.05. In addition, we imposed a sec- 
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scale factors of the Schechter function fit to the bright- 
and faint-end of the LFs, respectively. When iV-J" and N% 
are free to vary, we find that the geometric mean of the 
ratio (N$/N*) = 2.12 when measured over the four ra- 
dial bins from (r/r2oo) = 0.0 to 1.0. Due to the relatively 
bright absolute magnitude limits (—16.5) and the strong 
coupling of the Schechter function parameters, we have 
chosen to derive the faint-end slope a as the primary 
parameter of interest and thus set {N^/Nf) = 2. 

The resultant Schechter function fit parameters for the 
composite LFs are tabulated in Table 2, along with the 
1(7 uncertainties, where M£ and M£ are the turnover 
in the bright and faint Schechter functions, respectively. 
Since the faint-end slope of the first function has been 
fixed at a\ = —1.0, only a-i is presented. The individ- 
ual clusters that comprise the sample of 42 clusters used 
to construct the composite LF cover various fractions of 
(r/r2oo)- This results in a different number of clusters 
contributing to the LF counts in each of the four sepa- 
rate radial bins (see column 7 from Table 2). 

Examination of Figure [3] and Table 2 clearly indi- 
cates that the faint-end slope tends to become steeper 
as cluster-centric distance increases. The faint-end slope 
is significantly steeper, as measured out to (r /r^oo) = 1, 
than the traditional value of a = —1.25 (Schcchtcr 
1976). To facilitate the comparison between the dif- 
ferent LFs, the LF for each of the four radial bins has 
been plotted together in Figure [4] The outer three LFs, 
0.2 < (r/r2oo) < 1.0, have been scaled to match the 



counts of the inner-most LF in the —22 < Mr c < —21 
magnitude range. This figure clearly demonstrates the 
trend of a steepening of the faint-end slope with increas- 
ing cluster-centric radius. 

To determine what effect imposing the constraint 
(N£/N*) = 2 has on the robustness of the measured 
radial-dependent change in the faint-end slope, we refit 
our LFs by allowing N* and to vary. Due to the 
degenerate nature of fitting for A^*, M* and a simul- 
taneously, we measure the ratio (N^/N*) for the total 
cluster sample from four annuli by; i) fixing ai — —2.01, 
M* = —22.31 (mean values averaged over radial bins) 
and a.\ = —1.0, and then solving for ATf , A^ and M|, 
and ii) fixing M 2 * = -18.04, M{ = -22.31 (mean values 
averaged over radial bins) and ol\ = —1.0, and then solv- 
ing for AT-j* , AT| and ct2 ■ The increase in the steepness of 
the faint-end slope with increasing cluster-centric radius 
will be manifest by an increase in the (A^ /N*) ratio with 
increasing radius. 

For case (i) we find that the ratio (A^/A^*) = 
(0.49,0.50,5.10,13.00) for (r/r 200 ) = (0.0 - 0.2,0.2 - 
0.4,0.4 — 0.6,0.6 — 1.0). For case (ii) we measure 
(JVJ/ATf) = (0.85,2.76,2.60,4.19) for (r/r 200 ) = (0.0- 
0.2, 0.2 - 0.4, 0.4 - 0.6, 0.6 - 1.0). These results are con- 
sistent with the increase in the steepness of the faint-end 
slope with increasing radius. Due to the degeneracy in 
fitting the LF, the relative fraction of bright and faint 
cluster galaxies can be best studied by integrating over 
the respective LFs. This will be discussed in detail in 
the context of the dwarf-to-giant ratio of red- and blue- 
selected galaxies in Paper IV of this series. 

We compare our composit e LF with two recent stud- 
ies by converting to R c using Fukugita et a . (1995) and 



adopting our distance scale. IGonzalez et al 



((2006| used 

a double Schechter function to characterize the compos- 
ite LF of 728 groups/clusters selected from SDSS DR3. 
They find M* = -23.12 ± 0.12 with a = -1.89 ± 0.04 
for galaxies selected within 1.0 Mpc of the group cen- 
ter. Dividing their sample into two radial bins (0.0- 
0.6 and 0.6-1.2 Mpc), Gonzalez et al. determined that 
M* gets brighter with increasing radius (—22.9 ± 0.2 to 
—23.5 ± 0.2), while the faint-end slope becomes steeper 
{a = -1.80 ± 0.03 to a = -1.99 ± 0.03). The mea- 
sured range in M* for this study is significantly brighter 
(~ 0.6 — 1.1 mag) than our results tabulated in Table 
2. This may be related to the fact that the Gonzalez et 
al. sample consists predominately of group systems. For 
the faint-end, the measured slope values are consistent 
with our results, including the trend for a steepening of 
a with increasing radiu s . 

Pop esso et al.l (|2006h studied the composite LF of 69 
clusters from the ROSAT All Sky Survey/SDSS galaxy 
cluster catalog (RASS-SDSS) using the sum of two 
Schechter functions. The faint-end slope was found to 
increase from a = -2.02 ± 0.06 to -2.19 ± 0.09 when the 
sampling radius increased from r^oo to r2oo- This com- 
pares well with our measured range in a even though 
we use annuli to determine a rather than a circular re- 
gion with a changing radius. For the bright-end LF, 
Popesso et al. finds a varies from —1.05 ±0.07 for r5oo to 
— 1.09 ± 0.09 for r2oo- This result is consistent with our 
use of a fixed a\ = —1.0. Also, M* of the bright-end was 
found to brighten slightly from Mr c = —22.54 ± 0.13 
at ^500 to Mr c = —22.64 ± 0.16 for r2oo- These re- 
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suits are consistent within la of our measured values of 
-22.26±0.06 and -22.38±0.15 for the inner- and outer- 
most radial bin, respectively. In addition, we also detect 
a slight brightening of M* with increasing cluster-centric 
radius as reported by Popesso et al. for their sample. 

4.2.2. Red Sequence Luminosity Function 

To explore the dependency of the composite LF on the 
physical properties of cluster galaxies such as color, the 
cluster galaxy catalogs have been divided into various 
sub-samples that populate different regions of the CMD. 

The composite LF of galaxies located on the red se- 
quence of the CMR was constructed by selecting galax- 
ies using the following criteria: 1) galaxies are se- 
lected if they are within ±3.0 times the average dis- 
persi on of the Gaussian fit to the CMR (i.e., ±0.22 
mag; lLopez-Cruz et al.ll2004t ) or, depending on which is 
greater, 2) within 2.5 times the average uncertainty in 
the galaxy B — R c color redward of the red sequence fit 
(see Figure [5]). In practice, galaxies at the bright-end of 
the red sequence will be selected according to the first 
criterion; as fainter galaxies are chosen, the uncertainty 
in the B — R c color will increase to a value where galaxies 
will be selected via criterion 2. 

The first criterion defines the loci of the red sequence 
in the CMD. The second criterion helps insure that faint 
red sequence galaxies near the completeness limit are se- 
lected, even though they will be scattered further from 
the CMR than the limit imposed by criterion 1 due to the 
increased uncertainty in the measured color. As stated in 
the second criterion, only galaxies redward of the CMR 
were chosen. This was implemented to minimize the in- 
clusion of galaxies blueward of the red sequence (what 
we refer to as the "blue" galaxy population), which are 
scattered into the red sequence region due to the color 
uncertainty at faint magnitudes. To compensate for the 
"missing" members of the red galaxy population that 
are located blueward of the red sequence, the net galaxy 
counts for those galaxies fainter than the magnitude at 
which objects are culled based on criterion 2, have been 
increased by a factor of two, since in the absence of a blue 
galaxy population one might expect the red sequence 
members to be symmetrically distributed on cither side 
of the CMR. 

Figure [5] illustrates the region of the CMD where red 
sequence galaxies have been selected for Abell 260. For 
galaxies brighter than R c = 18.8 (indicated by the verti- 
cal dashed line), the 2.5<jb-r c is less than 3.0 times the 
average dispersion of the red sequence (0.22 mag), thus 
galaxies are chosen according to criterion 1 . For galaxies 
fainter than R c = 18.8, the uncertainty in the color mea- 
surement invokes selection by criterion 2. As depicted in 
Figure O only galaxies redder than the red sequence are 
selected for magnitudes fainter than R c = 18.8 mag. 

The net red sequence galaxy counts were calculated by 
subtracting the background field counts from the "raw" 
red galaxy counts. The background field galaxies were 
selected using identical color cuts as imposed on the red 
sequence galaxy sample. 

The composite red sequence galaxy LF for four differ- 
ent radial bins (same annuli as that used for the total 
composite LF) is depicted in Figure [6] The LFs have 
been constructed by scaling and combining clusters, com- 
plete to Mr c — —16.5, using the same method as that 



described for the total composite LF (§4.2.1). Visual 
inspection of the red sequence LFs shows that a single 
Schechter function, in general, would provide a reason- 
able fit to the overall shape of the LF except for the 
faintest magnitude bins. To compare the red sequence 
LF to the total LF, we fit the sum of two Schechter func- 
tions to the red sequence LFs. By allowing and iV^ 
to vary, we find that the geometric mean of the ratio 
(N£ /JVj")=0.4, averaged over the four radial bins. We 
thus impose that (iV|/JV*)=0.4 and fit for Mj*, M|, and 
a 2 ■ These fits are depicted in Figure [6] and the best-fit 
parameters tabulated in Table 3. We note that the LF 
fit to the faint component for the inner-most region is 
poorly constrained due to the lack of data points from 
the rising faint-end slope. We also find that when holding 
M 2 * fixed and allowing (AT| /ATf) and a 2 to vary, (JV| /N*) 
increases with increasing cluster-centric radius. 

In Figure [7] the four red sequence LFs have been su- 
perimposed by scaling the outer three LFs to match the 
inner-most LF in the magnitude range —22 < Mr c < 
—21. All four LFs appear to be equivalent to each other 
for magnitudes brighter than Mr c ~ —18. For mag- 
nitudes fainter than Mr c ~ —18, the inner- most LF 
appears to contain the smallest contribution from the 
faint galaxy component. In addition, a dip in the LF at 
Mr c ~ —18 is apparent for all four radial bins (see also 
Figure [6} . 

The red sequence LFs depicted in Figure [7] are remark- 
ably similar in shape to the composite LF for galaxies 
selected as having an r 1/,4 -like surface brightness profile 
given in lYagi et ail (|2002h (see their Figure 7) . This is 
not unexpected since the red sequence is dominated by 
early-type galaxies whose surface brightness distribution 
can be approximated by an r 1/,4 -like profile. The Yagi et 
al. composite LF is sampled to M Ra ~ —16.5, for our 
adopted cosmology, and thus compares directly with the 
red sequence LFs presented in this paper. An upturn in 
the red sequence LF at the faint-end is present in our 
sample (see Figure [6] and [7]) and in the composite r 1 / 4 
LF from Yagi et al. This upturn may be the result of 
scattering from the blue galaxy population; however, if 
this is the case, the scattered galaxies must also contam- 
inate the sample of r 1 / 4 -like surface brightness profiles 
measured by Ya gi et al., which se ems unlikely. 

In addition to lYagi et alJ (|2002f ). we also compare our 
red sequence LFs w ith the 2dF Galaxy Red shift Survey 
cluster sample from lDe Propris et al.l ((20031 ). This study 
presents composite LFs for a sample of 60 clusters based 
on the selection of galaxies relative to a classification pa- 
rameter derived from the principal com ponent analysis 
of spectra from iMadgwick et~afl (|2002t ). LFs are con- 
structed for early-, mid-, and late- type classes, with spec- 
tral classification based on star formation rates rather 
than morphological type. Assuming that the early-type 
spectral class is associated with red sequence cluster 
galaxies, we compare our results from Table 3 with De 
Propris et al. A single Schechter function fit to the early- 
type LF yields M* j = -20.04±0.09 and a = -1.05±0.04 
for Mb j < — 15. Converting to R c and using our distance 
scale, we find M* Rc = -22.92±0.09for M Rc < -18. Since 
De Propris et al. samples galaxies out to 3.0 Mpc, we 
compare our value for M* measured in the outer- most 
radial bin. From Table 3 we have M Rc = -22.61 ± 0.19, 
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which is within 1.5c of the De Propris et al. result. The 
De Propris et al. faint-end slope, a — —1.05 ± 0.04, is 
consistent with the value of a = —1.0 that we assume 
for the bright-end Schechter function given that De Pro- 
pris et al. only samples to Mr c ~ —18. Examination of 
our Figure [7] indicates that a becomes steeper for magni- 
tudes fainter than —18. Thus, De Propris et al. does not 
probe faint enough to sample the increasing faint galaxy 
component. 

In §4.2.1 we comp ared our total composite LF with 
iPopesso et al.l (|2006|) . That study also presents the com- 
posite LF for early-type cluster galaxies based on se- 
lecting galaxies with 2.22 < u — r < 3. Converting to 
R c and our distance scale, Popesso et al. finds M R = 
-22.27 ±0.14, with a faint-end slope of a = -2.01 ±0.11 
for galaxies measured within r2oo- Our M* values range 
from -22.28 ±0.07 to -22.61 ±0.19, thus we are in good 
agreement with Popesso et al. given that we measure 
LFs in a series of concentric annuli. The faint-end slope 
of Popesso et al. is flatter than our values tabulated in 
Table 3 (—5.26 < a < —2.83). This discrepancy is most- 
likely due to the poor fit at the faint-end of our LFs as 
indicated by the large reduced-x 2 values in Table 3. In 
fact, Popesso et al. samples > 1 mag deeper than our 
data set, and thus is able to place a better constraint 
on the faint-end slope. We also note that the rise in 
the faint-end slope of the early-type LF from Popesso et 
al. occurs at approximately the same magnitude as our 
study (M Rc 18.5). 

4.2.3. Blue Galaxy Luminosity Function 

We construct the LF of blue cluster galaxies by se- 
lecting galaxies in the CMD blueward of the red galaxy 
sample. The selection of the galaxies is illustrated in 
Figure [3 using Abell 260 as an example. Here, galaxies 
brighter than R c — 18.8 are designated as blue galaxies if 
their B — R c color is bluer than the blueward limit used 
to select the red sequence population. The color criterion 
for blue galaxies at M* is B - R c = 1.4 for Abell 260. 
This corresponds to an Sab-type galaxy usin g B — R c 
galaxy colors tabulated in lFukugita et all (|1995[ ). Galax- 
ies fainter than R c = 18.8 mag are selected if they are 
located in the region of the CMD that is bluer than the 
area used to select the red sequence galaxies and brighter 
than the completeness magnitude. An additional correc- 
tion to the galaxy counts is made by subtracting the net 
red sequence galaxy counts (measured from the region 
of the CMD fainter than R c — 18.8) from the net blue 
cluster galaxy counts. 6 This correction is necessary to 
account for the red sequence galaxies that inhabit the 
region in the CMD that is bluer than the area used to 
select the red galaxy population at the faint-end. 

The blue galaxy LF for four different annuli (same ra- 
dial bins as that used for the total and red composite 
LFs) is depicted in Figure [9] The LFs have been con- 
structed by scaling and combining clusters, complete to 
Mr c = —16.5, as described previously for the total and 
red sequence LFs. As determined for the total and red 
sequence LFs, a single Schechter function is unable to 

6 The net blue cluster galaxy counts were calculated by sub- 
tracting background field galaxies, selected using the same color 
criteria as the blue cluster galaxy population, from the raw blue 
cluster galaxy counts. 



adequately describe the shape of the blue galaxy LF. A 
sum of two Schechter functions (solid line in Figure [9]) is 
fit to the blue LFs, with the best-fit Schechter parame- 
ters given in Table 4. For the LF fit procedure we set 
(AVf /-/V-j*)=3, which is equal to the geometric mean of the 
ratio averaged over all four radial bins when jV-j" and A^jf 
are allowed to vary. 

In Figure [10] the blue LFs have been superimposed by 
scaling the outer three LFs to match the inner-most LF 
in the magnitude range —22 < Mr c < —21. In gen- 
eral, all four LFs appear to be very similar in shape, 
with a slight tendency for a steeper faint-end slope at a 
larger cluster-centric distance. This trend can be ascer- 
tained by examining the value of oti measured for the 
four LFs (see Table 4). We note that if we fix M| and 
solve for (N%/N?) and a 2 , (N%/Nf) decreases slightly 
with increasing cluster-centric radius. We also note that 
the blue LFs are similar in shape to the composite LF 
composed of galaxies h aving an exponen tial-like surface 
brightness profile from lYagi et a l. (2002) (see their Fig- 
ure 9). 

Anal ogous to the comparis on of our red seque n ce LF 
with iDe Propris et alJ (|2003h and IPopesso et all (|2006l ) 
in §4.2.2, we also compare our blue LF with the late- 
type LFs presented in these studies. Correcting for filter 
and cosmology difference, the single Schechter function 
fit to the late-type LF (Type 3+4) from De Propris et 
al. yields M* Ra = -22.02 ± 0.19 with a = -1.30 ± 0.10. 
This value for M* agrees well with our measured range 
(-21.96 < M* < -21.81), while our faint-end slope 
is much steeper (—1.62 < a < —1.82). As described 
in §4.2.2, the De Propris et al. study does not probe 
the faint-end of the LF to the same depth as our study 
(AM ~ 1.5 mag). The fit of the faint-end slope is fur- 
ther complicated by the large completeness corrections 
required for the faintest two magnitude bins. Thus, it is 
not surprising that the faint-end slope of the late-type 
LF from De Propris et al. is flatter than our value tab- 
ulated in Table 4. The trend of a steeper a for late-type 
galaxies compared to early-type systems, is consistent 
w ith our results for blu e versus red galaxies. 

IPopesso et al.l (|2006h presents the late-type composite 
LF constructed by selecting galaxies with u — r < 2.22 
for a sampling radius of r 2 oo- Unlike our results for the 
blue LF, a single Schechter function is found to be an 
adequate fit to the late-type LF. Converting the SDSS 
r-band to R c and employing our distance scale, the late- 
type LF from Popesso et al. has M* Ra = -23.41 ± 0.52 
and a = -1.87 ± 0.04. The value of M* from Popesso 
et al. is significantly brighter (~ 1.5 mag; 2.7cr level) 
than the average value measured for our blue LF. The 
Popesso et al. result is also ~ 1.4 mag brighter (2.5a 
level) then De Propris et al. for their late-type LF, and ~ 
0.8 mag brighter than the exponential composite LF from 
lYagi et all (|2002j ). We also note, based on De Propris et 
al., Yagi et al., and this study, that M* is brighter for 
the red (early-type, r 1 / 4 -like) LF than for the blue (late- 
type, exp-like) LF. This is not the case for Popesso et 
al., where the early- type LF has a fainter M* (~ 1.1 
mag) than the late-type LF. We do not understand this 
discrepancy, but speculate that it may be due to the 
selection of galaxy types by Popesso et al. based on u — r 
color. The faint-end slope of the Popesso et al. late-type 
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LF (a = —1.87 ± 0.04) is consistent with the range in 
slope values measured for our blue LF (—1.62 ± 0.05 to 
-1.82 ±0.10). 

In Figure [TT] the total, red sequence, and blue LFs for 
the four radial bins are directly compared. To aid this 
comparison, the red sequence and blue LFs have been 
matched to the total LF in each radial bin by scaling 
the red and blue galaxy counts such that the number 
of galaxies between —24 < Mr c < —16.5 is equal to 
the total LF in the same magnitude range. Inspection 
of Figure [TT] shows that the total and red sequence LFs 
are composed of a bright galaxy component that is much 
more significant than the contribution from the blue LFs. 
For the faint-end, the blue galaxy sample appears to con- 
tribute the greatest fraction of the faint dwarf galaxy 
population, with the largest difference apparent for the 
inner-most radial bin where (r/r2oo) < 0.2. 

It is also interesting to observe that M* for the faint LF 
component gets brighter with increasing cluster-centric 
radius for the red galaxy population and fainter for the 
blue galaxies. We will fully discuss this feature in the 
context of color-selected giant and dwarf galaxies in Pa- 
per IV of this series. 

4.3. Bright-End of the Luminosity Function 
4.3.1. The Distribution of M* 

The variation in the value of M* of the cluster LF can 
be used to gauge whether the LF is universal in shape at 
the bright-end, or whether luminosity segregation takes 
place (i.e., bright galaxies occupy preferentially the cen- 
tral regions of clusters). The measurement of M* has 
been important, hi storically, for potential use as a "stan- 
dard candle" fe.g. . lSchechQfl976l : lDresslerlll978l : iCollessI 
1989) and quantifying any variation in its value is impor- 
tant. We adopt the premise that the bright-end of the 
LF is universal if the median uncertainty of M* is com- 
parable to the measured dispersion of the distribution. 

To measure the variation in the bright-end of the LF, 
the value of M* was determined for our sample by fit- 
ting each cluster with a single Schechter function as de- 
scribed in §4.1; i.e., fitting the net galaxy counts in the 
magnitude range —24 < Mr c < —20, with the faint-end 
slope fixed at a = —1. Only galaxies measured within 
a cluster-centric radius of (r/r 2 oo) < 0.4 were selected. 
This allows us to maximize the number of clusters in 
our sample while normalizing each one to the same dy- 
namical radius without being adversely affected by small 
number statistics. Values for M* for 55 clusters satis- 
fying the magnitude and radius criteria are tabulated in 
Table 1. 

In Figure [12] we plot the distribution of MJj , which 
appears approximately Gaussian with a mean {M^ ) — 
—22.24 ± 0.06 and a dispersion of a = 0.31 mag. The 
median uncertainty in the measurement of for all 
55 clusters is 0.38 mag, thus supporting the universality 
of M * within the m e asure d range of AM* ~ 0.3 magni- 
tude. iPopesso et alj (|2005[ ) presents the histogram distri- 
bution of M* for LFs measured within a radius of 2h^Q 
Mpc. Their distribution is Gaussian-like in shape with 
an estimated mean of ~ —22.6 and a dispersion of ~ 0.5 
mag (transformed to our filter and distance scale). Given 
that the Popesso et al. LFs are measured with a fixed 
aperture, and at a larger radius than ours, it is not un- 



expected that they obtain a somewhat larger dispersion 
and a brighter mean than our results, which may 
arise from a mild dependence of M* on cluster-centric 
radius (see §4.3.2). 

Given the size of the dispersion, our results are not 
too dissimilar even though Popesso et al. used a fixed 
physical length scale for the counting aperture size rather 
than scaling relative to a dynamical radius like r2oo- 

Since the uncertainty in is expected to be a func- 
tion of cluster richness (assuming Poisson statistics), we 
measure the median uncertainty and dispersion of M* 
for three bins based on cluster richness; i) B gc > 1500, 

ii) 1000 < B gc < 1500, and iii) B gc < 1000. The di- 
viding B gc values are equivalent to line-of-sight veloc- 
ity dispersions of a v ~ 960 km s " 1 for B gc = 150 and 
a v ~ 750 km s- 1 for B gc = 1000 (|Yee fc Ellingsonll2003h . 

For clusters with B gc > 1500, we find a median un- 
certainty in MJj of 0.22 mag with a corresponding 
dispersion of 0.25 mag. For the cluster sample with 
1000 < B gc < 1500, we obtain a median uncertainty 
of 0.35 mag and a dispersion of 0.35 mag. Finally, for 
the poorest group we measure a median uncertainty of 
0.51 mag and a dispersion of 0.32 mag. A comparison 
of the median uncertainty and dispersion of for our 
three cluster sub-samples supports the hypothesis that 
the bright-end of the cluster LF is universal at the ~ 0.3 
mag level. 

A direct comparison of (M^ ) with other published 
values is problematic gi ven that is not independent 
of a (|Schechter|[l976l ). iPiranomonte et alj (|200lD finds 
M* = -22.02 ± 0.16 for the combined LF from a sample 
of 80 clusters using a Schechter function fit with a — 
— 1.01. The composite LF for this study was constructed 
by weighting clusters according to their richness. Trans- 
forming M* to the R c passband (jFukugita et al.l fl995h 
and our adopted cosmology, we find = —22.39 ± 

0.16, which is in very good agreement with our result. 
From a sample of 69 cluster s extracted from the RASS- 
SDSS, iPonesso et all (I2006T ) found M* = -20.94 ± 0.16 
and a = —1.09. The cluster galaxies were measured 
within a cluster-centric radius of r2oo- Transforming to 
our cosmology and passband yields = — 22.64±0.16, 
w hich is consistent wi th our measurement. 

lHansen et alj (|2005f ) presents M* for a sample of clus- 
ters/groups selected from the SDSS Early Data Release. 
The Hansen et al. study tabulates M* for several clus- 
ter samples of various richness with a faint-end slope 
fixed at a = —1.0 (see their Table 2). The value of 
M* for their richest sub-sample, which compares more 
directly with our Abell clusters, was determined to be 
M* = —20.86 ± 0.05. Transforming to our adopted cos- 
mology and i? c -band filter, we find = -22.53±0.05. 
This result is based on cluster galaxies measured within 
r2oo and is in excellent agreement with our results. 
Hansen et al. also suggests that M* is correlated with 
cluster richness in the sense that richer clusters have a 
brighter M*. For their three richest bins, M* differs 
by ~ 4cr between the poorest and richest sub-samples. 
To search for a similar correlation in our Abell sam- 
ple, we have divided our clusters into three bins based 
on Bgc] i) Bgc < 1000, ii) 1000 < B gc < 1500, and 

iii) Bgc > 1500. For these three sub-samples we find 
M* Rc = (-22.36 ±0.08, -22.25 ±0.07, -22.28 ±0.05) for 
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B gc = (< 1000, 1000 - 1500, > 1500). No correlation be- 
tween Afjj and cluster richness (as measured by B gc ) is 
evident. 

To search for a possible correlation with richness for 
individual clusters, we plot in Figure [13] the distribu- 
tion of Mfj c as a function of B gc . A S pearman rank- 
order correlation coefficient (|Press et al.lll992l) indicates 
that a correlation between MJj and B gc is not significant 
(r s = —0.08), with a 45% probability that the variables 
are correlated. It is not too surprising that cluster rich- 
ness, as measured by B gc , does not have a strong influ- 
ence on Mjj since we have normalized our cluster sample 
in terms of f2oo- In Figure [131 we do see a trend where 
the spread in the measured value of Mf, decreases with 
increasing B gc . This is most-likely due to the increased 
uncertainty in the measurement of for poorer clus- 
ters. For the ten clusters with the smallest value of B gc , 
the average uncertainty in MjjJ is 0.69 mag, while for 
the ten clusters with the largest value of B gc , the mean 
u ncertainty is .22 m ag . 

iHilton et al.l (|2005l ) finds a correlation between M£. 

and X-ray luminosity such that low-Lx clusters have a 
brighter M* than high-Lx systems (AM ~ 0.51 mag). 
To compare our results, w e use Lx measurements from 
lEbeling et~afl (|1996Ll2000T ) for a sub-sample of 41 clusters 
in common and construct composite LFs for a low-Lx 
(Lx < 3 x 10 44 ergs s _1 ; 22 clusters) and a high-Lx sam- 
ple (Lx > 3 x 10 44 ergs s _1 ; 19 clusters). Fitting a single 
Schechter function to the LFs yields M£ = -22.25±0.06 
for the low-Lx sample and M* Rc = -Tl.Tl ± 0.05 for the 
high-Lx group. Thus, there is no significant difference 
in M* when dividing our sample in terms of Lx- The 
discrepancy between our result and Hilton et al. may be 
due to the higher fraction of low-mass systems in their 
sample. For example, using the same dividing Lx as 
Hilton et al. (0.36 x 10 44 erg s _1 ), our low-Lx sample 
contains only 3 clusters (7% of the total), while the corre- 
sponding Hilton et al. sample contains 49 clusters (72% 
of the total) . 

To search for a correlation of L F parameters with 
cluster mass, iDe Propris et al.l (|2003ft divided their sam- 
ple into two groups based on velocity dispersion (a ^ 
800 km s _1 ). A Schechter function fit to the low- and 
high-er group yields the same M* at the ltr level. This 
is consistent with our findings that M* is not correlated 
with cluster mass. 

The characterization of clusters based on the rel- 
ative contrast of the brightest cluster galaxy defines 
the Bautz-Morgan type (BM-type) classification scheme 
(|Bautz fc Morgan! [1970). The evolutionary state of a 
cluster, as characterized by its BM-type, is expected to 
be correlated with M* if the BCGs are created via the 
merger of giant galaxies t hrough a cannibalism-like pro- 
cess (e.g.. [Dressier! I1978I ) . Galactic cannibalism would 
reduce the number of bright galaxies and therefore shift 
Mft to a fainter value. To test this scenario, we plot in 
Figure [TilM* versus BM-type for 54 clusters (A1569 has 
no published BM-type). The solid line in Figure [T3] de- 
picts a least-squares fit to the data and suggests that 
a weak correlation exits in the sense that early BM- 
type clusters have a fainter M* than later BM- types. A 
Spearman rank-order test gives a correlation coefficient 
of r s = —0.28, with a 96% probability that these two 



variables are correla ted. Thi s resul t lends support to the 
theoretical study of iMerritH (|1984l ) in which the forma- 
tion of the BCG occurs while the cluster is collapsing. 
It is expected that relatively little evolution of the BCG 
happens after the cluster is virialized. If we associate the 
early BM-type clusters with fully relaxed virialized sys- 
tems and late BM-types with unrelaxed systems in the 
process of collapsing, the trend of the correlation between 
M* and BM-type from Figure [14] can be explained. 

Since M* is correlated with BM-type, we correct the 
values of M* for this correlation and find that the disper- 
sion in Mft decreases from 0.31 mag to 0.24 mag. This 
result is consistent with the universality of M* within 
the 0.3 mag range. 

To investigate the possible correlation between M* and 
B gc for selected BM-types, we divide our 54 clusters into 
early BM-types (I and I II) and late BM-types (II, II-III, 
and III). A Spearman rank-order correlation analysis for 
the early BM-type sample yields r s — —0.27, with a 72% 
probability that M* and B gc are correlated. For the 
late BM-type clusters we find r s = —0.08 with only a 
39% probability that M* and B gc are correlated. These 
results indicate that no significant correlation exists be- 
tween M* and cluster richness among similar BM-type 
clusters when measuring galaxies within an equivalent 
dynamical radius and thus further supports the notion 
t hat M* is universal with in the measured uncertainty. 

IDe Propris et all ([20031 ) compares their composite LF 
fit parameters for two sub-samples divided according to 
BM-type. The difference in M* between the early BM- 
type group (I, II-II, II) and the late BM-type sample 
(II-III, III) is 0.06 ± 0.22 mag, where the uncertainties 
in M* are added in quadrature. The difference in M* 
between our BM-type I and III systems is AM = 0.3 
mag, which is significant at the 3<r level. Using identical 
BM-type bins as De Propris et al., we find AM = 0.2 
mag at the 3cr level. The discrepancy between our result 
and De Propris et al. may be related to differences in 
the technique used to construct the composite LFs. For 
example, De Propris et al. includes the BCG in their 
LFs. In addition, instead of using a dynamically scaled 
radius, De Propis et al. used all galaxies within a fixed 
aperture of 1.5h^ 00 Mpc, which may add to the dispersion 
of M*. Also, galaxy magnitudes from De Propris et al. 
are based on magnitudes, which are more susceptible 
to recent star formation and dust attenuation than R c 
magnitudes. In addition, De Propris et al. fits for the 
value of a while we impose the constraint that a = — 1. 
Since M* and a are correlated, a steeper a will in general 
yield a fainter M*. The faint-end slope of the early BM- 
type LF from De Propris et al. is flatter than the slope 
of the late BM-type LF. Forcing a to be the same will 
increase AM between these two samples. Although the 
difference in M* between the BM-type samples in De 
Propris et al. is not significant, it is interesting to note 
that the late BM-type sample has a slightly brighter M*, 
equivalent to the trend we measure. 

4.3.2. The Cluster- Centric Radial Dependence of M* 

The variation of M* as a function of cluster-centric 
radius was examined by measuring the bright-end of the 
composite LF constructed from our entire 57 cluster sam- 
ple. These clusters are 100% photometrically complete to 
Mr c — —20, and thus maximize the number of clusters 
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used to determine M* in order to minimize its measured 
uncertainty. (Recall that the M* values tabulated in Ta- 
bles 2-4 are measured for a sample of 42 clusters that 
are complete to Mr c — —16.5.) 

The dependence of M* on cluster-centric radius was 
determined by fitting a single Schechter function to the 
net galaxy counts from —24 < Mr c < —20 for the four 
radial bins used previously (for example, see Figure [3]). 
The faint-end slope of the Schechter function was fixed at 
a = —1.0 as was done when measuring M* for the indi- 
vidual cluster LFs. Since different clusters cover various 
fractions of r2oo, the total number of clusters contribut- 
ing to the composite sample will vary with cluster-centric 
radius. 

In Figure [15] we plot M* as a function of (r/r2oo) 
for the total composite cluster sample, as well as sam- 
ples compiled by selecting red and blue galaxies from 
our 57-cluster sample (see §4.2 for color selection crite- 
ria). Inspection of Figure [T5] reveals that for the red 
and blue composite samples, M* gets brighter with in- 
creasing cluster-centric radius. For the total composite 
sample, we find a similar correlation, although the trend 
is much weaker. This weaker trend is due to the in- 
terplay between the relative dominance of the red and 
blue galaxies as a function of cluster-centric radius. M* 
near the center is dominated by red galaxies; whereas at 
larger radii, blue galaxies with a fainter Mr c become an 
increasingly more important component, flattening the 
M^ dependence on radius. The difference in M* be- 
tween the inner- and outer-most radial bin is AM — 0.48 
mag for the red sequence LF (3. 3c level) and AM — 0.27 
mag for the blue LF (1.2cr level). 

To quantify the correlation for each of our three com- 
posite cluster samples, we calculate the l inear correlation 
coefficient r statistic ([Press et al.lll992t ). For the total 
composite cluster sample, we find that the linear corre- 
lation coefficient is r = —0.91, with a 92% probability 
that M* and radius are correlated. The red sequence 
composite cluster sample yields r = —0.98, with a 98% 
probability of a correlation. For the blue galaxy cluster 
sample, the linear correlation coefficient is r = —0.98 and 
a 98% probability of a correlation. 

5. DISCUSSION AND CONCLUSIONS 

5.1. Universality of the Luminosity Function 

We have examined the distribution of M* for the 
bright-end of the LF in the core of Abell clusters. We 
found that the dispersion in M* is comparable to the av- 
erage measured uncertainty, even when dividing the clus- 
ter sample into different richness groups based on B gc . 
We use this result to indicate that M* is universal at 
the 0.3 mag level for a restricted magnitude range and 
when measured within a specific dynamical radius. In 
addition, we find a weak trend in which early BM-type 
clusters have a fainter M* than late BM-types. 

One of the primary goals of this paper is to explore 
the change in the cluster i? c -band LF as a function of 
cluster-centric radius. From a sample of 57 low-z Abell 
clusters, we have measured LFs covering various frac- 
tions of r2oo for both the red sequence and blue cluster 
galaxy populations. Our results indicate that the over- 
all shape of the LF is dependent upon distance from the 
cluster center. In general, the LFs exhibit an increase in 
the steepness of the faint-end slope with increasing ra- 



dius. The radial dependence of the rate-of-change in a 
is greatest for the total cluster sample (i.e., red plus blue 
galaxies), while the blue LF is less dependent. The red 
sequence LF is mostly flat (a ~ —1) over a span of ^ 5 
magnitudes (—23 < Mr c < —18), with a rising faint-end 
for Mr c > —18. In contrast, the blue galaxy LF is much 
steeper (a ~ —1.8) over this same luminosity range, with 
minimal change in shape out to ^200- The very rapid 
increase in a for the total LF is likely due to a com- 
bination of steepening slope for both red and blue LFs 
and the increasing dominance of the blue population at 
large cluster-centric radii. In addition, the red sequence 
LF has a much brighter characteristic magnitude (~ 0.6 
mag) over all radii than the blue galaxy luminosity dis- 
tribution. These results lend support to several recent 
studies that have observed similar characteristics for red 
and blue galaxy LFs drawn from low- to high-density en- 
vironments (e.g., iBeiiersbergen et al.ll2002at iGoto et ail 
2002 5 IBaldrv et al.H2004fl . 

The general trend for the LF to become flatter with 
decreasing cluster-centric radius supports the hypothesis 
that dwarf galaxies are tidally disrupt ed near the cluster 
center . This idea has been used by lLopez-Cruz et all 
(|1997f ) to help explain the formation of BCG halos and 
the origin of a large fraction of the gas content in the 
intracluster medium. 

The dependence of the shape of the LF on cluster- 
centric radius provides strong evidence that the relative 
mixture of giant and dwarf galaxies depends on the frac- 
tion of the virial radius that is measured. This argues 
against the global universality of the cluster LF for the 
magnitude interval —26 < Mr c < —16.5 and suggests 
an environmental influence. The non-parametric galaxy 
dwarf-to-giant ratio will be explored in Paper IV of this 
series. 

In FigureHTOlwe plot composite cluster LFs from several 
published sources 
Goto et all 12001 



Collessll989tlPiranomonte et al.ll2001l: 
De Propris et al.l 120031 : lHansen et al.1 



20051: iPopesso et al.l 120 061) . These LFs h ave b een trans 
formed to Mr^ using iFukugita et"aH (fl995h and our 
adopted cosmology. Also depicted in Figure [TH] are the 
total composite LFs from this study for the inner- and 
outer-m ost radial bin. In add ition, the SDSS field galaxy 
LF from Blan ton et al.l (|2003f ) is included for comparison 
purposes. 

Figure [16] illustrates that the slope of the LFs are very 
similar at bright magnitudes where the giant galaxies 
have the greatest influence on the shape of the LF. The 
main difference in the slope arises at the faint-end where 
the influence of the dwarf galaxies tends to increase the 
steepness of a. The sampling depth and effective cluster- 
centric radius thus has a major influence on the measured 
shape of the cluster LF since the inclusion of different 
amounts of the dwarf galaxy population will directly im- 
pact the faint-end slope. Thus, the evidence supports the 
notion that the cluster LF is not universal in shape in a 
global sense. In addition, we measure a dispersion of 0.3 
mag in the value of M* for the depicted cluster LFs; the 
faintest value is Mjj = —22.26 for the inner-most annuli 
from this study and the brightest is Mjj = —23.14 from 
IGoto et alJ (|200l . 

5.2. Cluster Galaxy Population Gradients 
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As discussed in §4.3.2, Figure [TBI shows that the bright- 
end M* becomes brighter with increasing cluster-centric 
radius for both the red, and with a lesser significance, 
the blue cluster galaxies. The observed dimming of M* 
toward the cluster center can be explained as a simple 
fading of the galaxy population. For the blue cluster 
galaxies, the truncation of star formatio n as field galax- 
ies fa l l into the cluster environment fe.g.. lAbraham et al.l 
[19961: lEllingsonl l2003h would lead to a fading of M* 
with decreasing cluster-centric radius. Since clusters 
are believed to have formed via the i nfall of galaxie s 
along filamentary-like structures (e.g., iDubinskil fl998L 
it is expected that a population of infalling field galaxies 
can be detected in addition to the older, mainly early- 
type red galaxies. If star formation for infalling spi- 
ral ga laxies, which domina te the field galaxy population 
(e.g., iBinggeli et al.l fl988l ). is truncated via some type 
of dynamical mechanism (e.g., ram pressu r e stripping , 
galaxy harassment, e tc.; iMoore et al.l 119981 : iGotol [20051 : 
iRoediger et~aH 120061 ) . then a roughly continuous infall 
(but with an allowed variable rate) would be observed 
as a fading of the population toward the cluster center. 
Given enough time, infalling spiral galaxies may acquire 
characteristics that are similar to SO galaxies. In fact, 
this type of mechanism for SO formation has been pro- 
posed by numerous authors (see lDressler et al.lll999L and 
references therein), although the formation of field SOs 
has remained problematic for these models. 

The dimming of M|, or the decrease of (N%/Nf) 
for the faint-end of the red sequence LF for decreas- 
ing cluster-centric radii, places constraints on the evolu- 
tionary path of the faint blue galaxies. If these galax- 
ies simply fade and turn red, we would expect them 
to contribute to the red sequence LF by increasing the 
number of faint red galaxies in the central cluster re- 
gion. We thus suggest that the faint blue cluster galax- 
ies are destroyed in the central cluster region. This is 
not a far-fetched hypothesis since the faint blue galaxies 
are very similar to the low-mass dwarf spheroidal galax- 
ies, which are expected to unde rgo tidal disruption in 
the central cluster environme nt (Thompson fe Gregory! 
[19931: iGallagher fe Wvsdll99l iHilker et al.ll2003l) . These 
low luminosity blue galaxies could also be the source 
of the dwarf galaxies that get tidally disrupted and 
form the halo of BCG s as proposed in a model by 
lLopez-Cruz et all (|1997l ). We will explore this possibility 
further in Paper IV of this series. 

The fading of the bright-end M* for the red sequence 
galaxy population may be the result of two separate 
processes; the continuous fading of infalling red galax- 
ies (which may have turned red relatively recently due 
to the truncation of star formation) and galactic can- 
nibalism in the high-density central cluster region. For 
the infall scenario, red galaxies are expected to origi- 
nate from a population that had arrived in the cluster 
environment early in its history. The observed disper- 
sion of the CMR for the cluster red sequence places con- 
straints both on the formation timescal e of the early- 
type g alaxy population (i .e., z > 2: e.g.. IStanford et al.l 
19981 iDePropris et all [200l iGladders et alJ Tl999F 
Holden et all 120041 : lLopez-Cruz et all l200# and any 
episo des of recent star formation (e.g., iBower et alJ 
1998). Red galaxies which have entered the cluster envi- 
ronment during the earliest part of the cluster lifetime, 



would be expected to be the faintest since the time be- 
tween the last phase of major star formation and to- 
day would be the greatest for these galaxies (after the 
truncation of star formation, galaxies would evolve pas- 
sively with an associated reddening and dimming of their 
starlight). Under this scenario, we would expect that 
the blue galaxy population would show a greater rate-of- 
change in M* with radius due to the more recent decline 
in the star formation rate. The red galaxies would be 
expected to exhibit a more gradual change in M*, as 
compared to the blue galaxy population, since the red 
galaxies are just a passively evolving old galaxy popula- 
tion. 

As depicted in Figure [T5] the bright-end M* for the 
red sequence galaxies fades more rapidly toward the clus- 
ter center than for the blue galaxy population. Since a 
simple infall scenario as described predicts that the blue 
galaxy population should exhibit the more rapid decline 
in M* , we hypothesize that an additional mechanism 
may affect the rate-of-change of M* with radius for the 
red galaxy population. A clue to this additional effect 
is gleamed from Figure [T4l where the correlation between 
M* and BM-type is depicted for the total cluster sample. 
As described in §4.3.1, the trend for a brightening of M* 
with later BM-type can be explained by galactic canni - 
bali sm. As first theorized bylOstriker fe Hausmanl (|1977l ) 
and lHausman fc Ostrikerl (|1978f ) , as bright galaxies fall 
into the cluster center they will be swallowed by the gi- 
ant central galaxy, thus resulting in the fading of M* 
as bright galaxies are depleted from the LF. Since early 
BM-type clusters have, by definition, a bright central 
dominant galaxy, the correlation depicted in Figure [15] 
is expected. Galactic cannibalism will also result in the 
fading of M* with decreasing cluster-centric radius, as il- 
lustrated in Figure[T51 We performed a simple test of this 
hypothesis by calculating M * for our red cluster galaxy 
sample for early BM-type clusters (I and I-II) and late 
BM-type clusters (II, II-III, and III). We assume that the 
effect of galactic cannibalism on the radial dependence of 
M* will be greatest for the early BM-type sample. Com- 
paring M* between the inner- and outer-most radial bin 
shows that for the early BM-type sample, AM = 0.71 
mag with the two measurements different at the 2a level. 
For the late BM-type sample, we find AM — 0.36 mag at 
the 2ct level. Although the significance of the difference in 
M* for a given BM-type group is not high (partly due to 
the fact that only four clusters comprise the outer-most 
radial bin for the early BM-type clusters), the larger dif- 
ference in M* for the early BM-type sample suggests that 
the infall and galactic cannibalism scenarios may help to 
explain the observed rate-of-change of M* with radius 
for the red sequence galaxy population. 

In Paper IV of this series we will further elucidate 
the nature of the cluster galaxy population by exam- 
ining the radial-dependence of the dwarf-to-giant ratio 
and the blue-to-red galaxy count ratio. These results will 
complement the observations presented in this paper and 
provide additional constraints on the composition of the 
galaxy population in low-redshift Abell clusters. 
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APPENDIX 

A. QUANTIFYING SELECTION EFFECTS AND BIASES 

A major challenge when conducting any observational study is the ability to quantify the impact of selection effects 
and biases on the robustness of the results. In this appendix we explore several of these potential sources of systematic 
errors and quantify their affect on the derived galaxy cluster LF. 

A.l. Chance Projections 

A main concern in interpreting the results of a study based on the LF constructed from the statistical subtraction of 
background galaxies is the possible influence of projection effects. The chan ce projection of field galaxies, or unrelated 
groups, can appear as clusters when redshift information is not available (jValotto et al.l l200fl) . Selection of cluster 
members using color information (see §3.2) helps to alleviate some of this concern. To determine what impact chance 
projections may have on our LFs, we have div ided our cluster sample into two groups according to the number of 
published redshifts (e.g.. IStruble fc Roodlll999T ). The first group contains only those clusters which have at least 25 
spectroscopically-confirmed cluster members. The second group contains clusters in which the number of redshift- 
confirmed members is < 10. In this comparison it is assumed that clusters with a small number of confirmed members 
may still suffer from unknown projection effects. 

In Figure [17] we present LFs for the total composite cluster sample for four different radial bins. The cluster sample 
has been divided into two groups according to the number of confirmed cluster members based on the number of 
redshifts (i.e., 20 clusters with < 10 galaxy redshifts and 32 clusters with > 25 galaxy redshifts). The LFs for each 
radial bin have been scaled to match the LF for that particular cluster-centric radius from the total composite LFs 
depicted in Figure [3] The overall shape of the LFs generated for the two redshift-selected cluster groups are very 
similar in shape for all four annuli. Hence, this test indicates that the chance projection of background clusters/galaxy 
groups does not have a significant effect on our results. 

A. 2. Deprojecting the Luminosity Function 

The LFs presented in this paper are derived from the two-dimensional projected distribution of cluster galaxies. 
Since the full three-dimensional galaxy spatial positions cannot be resolved, the presence of contaminating galaxies 
from within the cluster can adversely affect the accuracy in determining the shape of the central cluster LF. Note that 
in this context, the projected galaxies are those which lie in the outskirts of the cluster and are projected onto the 
central region. This is in addition to the presence of fore/background galaxies that are unrelated to the target cluster, 
and are corrected using background corrections. The projection of galaxies located in the outskirts of clusters onto 
the central region will have its greatest impact on the faint-end of the central LF. This is due to the increase in the 
relative fraction of dwarf galaxies with increasing cluster-centric radius, as implied by the steepening of the faint-end 
slope with radius (see Figure EH) . 

Simulations by IValotto et al.l (120011 ) indicate that projection effects can severely affect the shape of the LF by 
artificially producing a steeper faint-end slope. Some of this steepening is due to the fact that cluster galaxies were 
not selected in terms of color (i.e., using the CMR), thus resulting in contamination by projected field galaxies. A 
certain portion can also be attribute d to the inclusion of galaxies in cluster outskirts that are projected onto the central 
region. iBeiiersbergen et"aT1 (|2002af ) published a study of the LF of the Coma cluster in which they corrected the LF 
for projection effects by subtracting the contribution of the outer Coma LF that is projected onto the central region. 
The resulting "deprojected" LF was measured to be marginally flatter than the projected 2-d LF. 

To determine the extent to which our LFs are affected by projection effects from the galaxy population in the cluster 
outskirts, we have measured the deprojected composite LF for the total, red sequence, and blue galaxy samples. We 
assume that the cluster galaxies are distributed symmetrically in a sphere about the cluster center with the method 
of deprojecting the LF depicted in Figure [TBI As shown, any sight-line to the central area of the cluster (region A) 
must pass through regions Bl and B2. Thus, the observed central LF will include galaxies that are located in regions 
Bl and B2. To correct for this effect, the LF determined from region C is subtracted from the LF calculated for the 
regions B1+B2+A. Before subtracting the two LFs, the LF measured for region C must be normalized to the same 
volume as that included in regions B1+B2. This "correction" factor is given by (Vi — V^)/(V3 — V%), where V\ is the 
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combined volume of regions B1+B2+A, V 2 is the volume of the inner region A (V 2 = (4.ir/3)Rf), and V 3 is the volume 
of the outer-most sphere (V3 = (47r/3)Ji|). The volume V\ can be calculated from V\ = (47r/3) [R\ - - i?^) 3/2 ] , 
where R\ is the radius of the inner sphere, and R2 is the radius of the outer sphere. 

Using this procedure, the deprojected composite LF for the total, red sequence, and blue galaxy samples were 
constructed for two radial bins, (r/r 2 oo) < 0.2 and 0.2 < (r/r 200 ) < 0.4. The large uncertainty in the net galaxy counts 
from the deprojection process limits our analysis of the LF to these two inner-most radial bins. In Figures [T51-I2T1 
comparisons between the projected and deprojected LFs for the three cluster galaxy samples are presented. For each 
figure the LFs have been scaled to match in the —22 < Mr c < —21 magnitude range to facilitate the comparison. 
The top panel for each of the three figures depicts the deprojected LF for the two inner-most radial bins. For all 
three galaxy samples, the LF for the (r/r2oo) < 0.2 radial bin (open triangles) has a flatter faint-end slope than for 
the 0.2 < (r/r2oo) < 0.4 radial bin (solid triangles). Thus, the trend of a steepening faint-end slope with increasing 
cluster-centric radius is valid for the deprojected LFs for the total, red sequence, and blue galaxy samples. 

In the middle panel of Figures [lT)H2"Tl the deprojected and projected LFs are compared for the (r/r2oo) < 0.2 
radial bin. For all three galaxy samples, the deprojected LF has a flatter faint-end slope than the projected LF. This 
demonstrates that galaxies from the cluster outskirts that are projected onto the central cluster region will result in 
an artificial steepening of the faint-end of the central LF. 

The bottom panel of Figures [T9M211 presents the deprojected and projected LFs for the 0.2 < (r/r 2 oo) < 0.4 radial 
bin. In general, the shape of the LFs are very similar given the size of the uncertainties in the net galaxy counts. No 
statistically significant difference in the shape of the deprojected and projected LFs can be discerned for either the 
total, red sequence, or blue galaxy cluster samples. 

The deprojection of the cluster LF indicates that the presence of galaxies from the outer cluster region can affect the 
slope of the faint-end, especially for the central region. The basic trend of a steepening faint-end slope with increasing 
cluster-centric radius is still evident from the deprojected LF data. 

A. 3. Faint Source Correction 

The shape of the faint-end of the cluster LF can be affected by bias as a result of counting ga laxies in the faint est 
magnitude bin. The faint galaxy correction (also known as the Eddington bias or correction; Eddingtonl Il940f ) is 
due to the fact that each observed galaxy has an associated magnitude uncertainty, causing galaxies to be scattered 
below and above our observed magnitude limit. If the number distribution of galaxies is identical over all magnitudes, 
the number of galaxies scattered above and below our magnitude threshold will be statistically equal. A cluster LF 
having a faint-end slope steeper than a = — 1 , will have a net number of galaxies scattered brighter than the observed 
magnitude limit. This will artificially enhance the number of detected galaxies at the faint-end of the LF. Since our 
galaxy counts are binned to produce the LF, the magnitude bin size relative to the uncertainty of the measured galaxy 
magnitudes will have a direct impact on the importance of this bias. 

To investigate the impact of the faint source correction on our measured LFs, simulated LFs using the shape 
parameters tabulated in Table 2 for the total composite LF were constructed. For each galaxy magnitude we include 
an offset calculated from the observed distribution of magnitude uncertainty aji c as a function of -Mr c , assuming 
Gaussian statistics. The magnitude offset applied to each simulated galaxy was randomly chosen from a Gaussian 
distribution with a dispersion of 1ctr c . Simulated LFs were assembled using the same selection criteria as that 
used for the observed LFs. The fractional change in the number of galaxies detected in the faintest magnitude bin 
(—17.0 < Mr c < —16.5) depends directly on the value of a 2 . Comparing simulated LFs to those tabulated in Table 2 
yields an increase from 0.0% to +0.4% in the number of detected galaxies in the faintest magnitude bin as the value of 
a 2 changed from —1.18 to —2.43. This demonstrates that the increase in the number of detected galaxies due to the 
scattering of faint galaxies into the faintest magnitude bin, is insignificant given the measured range in a 2 . This result 
is reasonable given the relatively large width of our magnitude bins (0.5 mag) and the average magnitude uncertainty 
(^0.1 mag) at our imposed faint-end magnitude limit. 

A. 4. Color Selection Bias 

The study of the LF of galaxies selected according to their position in the color-magnitude plane could be affected 
by a color bias. This bias is a result of scattering in B — R c and R c of galaxies due to photometric uncertainties. To 
understand the extent of this effect, an artificial galaxy catalog for Abell 260 was constructed by taking the original 
galaxy catalog and adding a small, random magnitude offset based on the measured R c magnitude and its uncertainty 
for each galaxy, assuming Gaussian statistics. By tracking the relative offset in position on the CMD between the 
original cluster galaxies and their simulated counterparts, the impact of the color bias can be measured. 

In Figure [22] a vector-type CMD for Abell 260 is presented for galaxies brighter than the cluster completeness limit 
(i? c = 20.7). The vectors trace the scattering path of a galaxy from its initial position in the observed color-magnitude 
plane to its position in the artificial catalog. For bright galaxies, the change in position is minimal compared to the 
size of the region from where the galaxies are selected (see Figures [5] and [8] for selection of red sequence and blue 
galaxies). The largest displacement in the color-magnitude plane occurs for the faint red galaxies where the relatively 
large uncertainties in B — R c produce a larger change in color. The magnitude of this displacement is not significant 
compared to the size of the regions used to select red and blue galaxies. This exercise justifies the method used to 
select faint red sequence galaxies (see §4.2.2, criterion 2). By restricting the selection of faint red sequence galaxies 
from a region redward of the CMR (see Figure 0), the contamination from blue galaxies scattered into the red sequence 
region is minimized. 
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A. 5. Dynamical Radius Variation 

As stated in Section 3.4, the value of r2oo determined from B gc (via equation 7) has an associated rms scatter of 
~ 15%. To determine whether a 15% scatter in r2oo will have a significant influence on our conclusions, a simulated 
composite LF was derived by randomly changing r 2 oo by ±15%. A comparison of the observed LF with the simulated 
LF is presented in Figure [33J for the four cluster-centric radial bins used previously. Inspection of Figure [53] shows that 
the two LFs are equivalent for each radial bin depicted. Thus, the 15% scatter in the value of r2oo is expected to have 
minimal impact on our results. This is at least in part due to the expected gradual change in the properties of the LF 
as a function of r2oo- 

A. 6. Exclusion of the Brightest Cluster Galaxy 

In §4.1 it was noted that the BCG was not included in the construction of the cluster LF since they do not 
appear, statistically, to be a natural extension of the LF. In Figure [24] the composite LF for the inner-most radial 
bin, (r/r2oo) < 0.2, is depicted for clusters (complete to Mr c — —20) with and without the inclusion of the BCG. 
As this figure demonstrates, the BCGs are not a simple extension of th e Schechter fun ction, which may indicate that 
BCGs are formed by a different process (e.g., mergers, cannibalism, etc.; Dressier [l 9 78t ) than the fainter cluster galaxy 
population. 
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CLUSTER PROPERTIES 
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A2244 


0.0997 


-22.14 


0.22 


1674 


243 


2.802 


0.420 


-17.5 


A2247 


0.0385 


-23.07 


0.78 


639 


151 


1.765 


0.265 


-16.5 


A2255 


0.0800 


-22.60 


0.21 


2278 


280 


3.248 


0.487 


-17.0 


A2256 


0.0601 


-22.51 


0.21 


2187 


274 


3.185 


0.478 


-16.5 


A2271 


0.0568 


-21.46 


0.51 


669 


157 


1.804 


0.270 


-16.5 


A2328 


0.1470 


-22.01 


0.20 


1935 


263 


3.004 


0.450 


-18.5 


A2356 


0.1161 


-22.25 


0.36 


964 


189 


2.150 


0.322 


-18.0 


A2384 


0.0943 


-22.38 


0.29 


1514 


232 


2.670 


0.400 


-17.5 


A2399 


0.0587 


-22.16 


0.54 


676 


157 


1.813 


0.272 


-16.5 


A2410 


0.0806 


-22.24 


0.71 


546 


145 


1.636 


0.245 


-17.0 


A2415 


0.0597 


-21.74 


0.44 


940 


184 


2.123 


0.318 


-16.5 


A2420 


0.0838 


-21.94 


0.27 


1239 


210 


2.425 


0.364 


-17.5 


A2440 


0.0904 


-22.34 


0.35 


1050 


196 


2.240 


0.336 


-17.5 


A2554 


0.1108 


-22.61 


0.35 


1221 


211 


2.408 


0.361 


-18.0 


A2556 


0.0865 


-22.47 


0.49 


796 


172 


1.961 


0.294 


-17.0 


A2593 


0.0421 


-22.25 


0.40 


1133 


200 


2.323 


0.348 


-16.0 


A2597 


0.0825 


-21.56 


0.58 


696 


163 


1.839 


0.276 


-17.0 


A2626 


0.0573 


-22.65 


0.50 


911 


181 


2.092 


0.314 


-16.5 


A2634 


0.0310 


-22.30 


0.32 


1109 


197 


2.299 


0.345 


-16.5 


A2657 


0.0414 


-22.44 


0.51 


723 


162 


1.872 


0.281 


-16.0 


A2670 


0.0761 


-22.52 


0.26 


1783 


249 


2.888 


0.433 


-17.0 



a The absolute R c magnitude represents our adopted 100% completeness limit. 
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TABLE 2 

COMPOSITE LUMINOSITY FUNCTION 



Radius 


Ml 


xl 


M 2 * 


«2 


xl 


No. of 


< Bg C > 


(r/r 2 oo) 


(Rc mag) 




(R c mag) 






Clusters 


(Mpc 1 - 8 ) 


0.0-0.2 


-22.26 ±0.06 a 


0.64 


-17.43 ± 0.07 b 


-1.45 ±0.10 


0.94 


42 


1066 ± 397° 


0.2-0.4 


-22.26 ±0.07 


1.04 


-18.22 ±0.05 


-1.81 ±0.04 


4.50 


39 


1052 ± 378 


0.4-0.6 


-22.36 ±0.10 


0.58 


-18.14 ±0.05 


-2.32 ± 0.05 


3.24 


28 


961 ± 347 


0.6-1.0 


-22.38 ±0.15 


0.67 


-18.39 ±0.06 


-2.46 ± 0.05 


2.99 


11 


683 ± 205 



a M£ is derived from a Schechter function fit to the bright-end of the LF with cti = —1.0. 



b M| and «2 are obtained from a Schechter function fit to the faint-end of the composite LF. 
c The uncertainty is calculated from the dispersion of the mean. 



Luminosity Function of Abell Clusters 



TABLE 3 

COMPOSITE RED SEQUENCE LUMINOSITY FUNCTION 



Radius 


Ml 


xl 


M* 


Ot 2 


xl 


No. of 


< Bg C > 


0/r 2 oo) 


(Rc mag) 




(R c mag) 






Clusters 


(Mpc 1 - 8 ) 


0.0-0.2 


-22.28 ±0.07 a 


0.69 


-16.95 ± 0.58 b 


-5.26 ± 15.51 


3.96 


42 


1066 ± 397 c 


0.2-0.4 


-22.36 ± 0.08 


1.29 


-17.81 ±0.22 


-3.30 ± 0.64 


1.45 


39 


1052 ± 378 


0.4-0.6 


-22.50 ± 0.12 


0.72 


-18.18 ±0.24 


-3.16 ±0.51 


1.16 


28 


961 ± 347 


0.6-1.0 


-22.61 ± 0.19 


1.14 


-18.60 ±0.36 


-2.83 ±0.53 


1.06 


11 


683 ± 205 



a Mj is derived from a Schechter function fit to the bright-end of the LF with oi = —1.0. 

b M| and 02 are obtained from a Schechter function fit to the faint-end of the red sequence LF. 

c The uncertainty is calculated from the dispersion of the mean. 
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TABLE 4 

COMPOSITE BLUE LUMINOSITY FUNCTION 



Radius 


Ml 


xl 


M 2 * 


»2 


xl 


No. of 


< Bg C > 


(r/r 2 oo) 


(Rc mag) 




(R c mag) 






Clusters 


(Mpc 1 - 8 ) 


0.0-0.2 


-21.96 ±0.26 a 


0.28 


-19.30 ± 0.10 b 


-1.62 ±0.05 


0.61 


42 


1066 ± 397° 


0.2-0.4 


-21.84 ±0.17 


0.36 


-19.28 ±0.09 


-1.64 ±0.05 


0.93 


39 


1052 ± 378 


0.4-0.6 


-21.81 ±0.21 


1.21 


-19.01 ±0.12 


-1.69 ±0.09 


1.05 


28 


961 ± 347 


0.6-1.0 


-21.87 ±0.27 


0.74 


-18.79 ±0.12 


-1.82 ±0.10 


0.61 


11 


683 ± 205 



a M£ is derived from a Schechter function fit to the bright-end of the LF with cti = —1.0. 



b M| and c*2 are obtained from a Schechter function fit to the faint-end of the blue galaxy LF. 
c The uncertainty is calculated from the dispersion of the mean. 
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Fig. 1a. — Solid points depict the i? c -band luminosity function for 57 Abell clusters. The solid line represents the best-fit Schechter 
function with a fixed faint-end slope of a = —1. The galaxy counts arc measured from within a cluster-centric radius of (r/r20o) = 0.4, 
except for A496 and A1142 due to the lack of adequate radial coverage. The brightest cluster galaxy has been omitted from each LF. 
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Continued 
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Fig. lc. — Continued 
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LogB gc (Mpc 18 ) 

Fig. 2. — Logarithmic correlation between r200 and B gc for 15 clusters from the CNOC1 cluster sample. The solid line is a fit using the 
BCES estimator to the 14 clusters depicted by the solid circles. The open circle represents the outlier cluster MS 1455+22, which was not 
used in the fitting process. The rms scatter in the derived values of r2oo is on the order of 15%. 



26 



Barkhouse et al. 




-25 -23 -21 -19 -17 -25 -23 -21 -19 -17 

Absolute R c Magnitude 



Fig. 3. — Composite total R c -band luminosity function for four cluster-centric annuli: A) (r/r2oo) < 0-2, B) 0.2 < (r/r2oo) < 0.4, C) 
0.4 < (r/r2oo) < 0.6, and D) 0.6 < (V/f"20o) < 1-0. The short dashed line represents a Schcchtcr function fit to bright-end with a fixed 
faint-end slope, a = —1. The solid circles depict the combined net galaxy counts for all contributing clusters. The long dashed line is a 
Schcchtcr function fit to the faint-end, while the solid line is the sum of the two Schechter functions. 
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Fig. 4. — Superposition of the total composite LF measured for four cluster-centric radial bins. The outer three LFs have been scaled 
to match the inner-most LF in the —22 < Mr c < —21 magnitude range. The plot symbols depict the following: solid circles — the net 
galaxy counts in the (r/r2oo) < 0.2 annulus; open squares — the counts in the 0.2 < (r/r2oo) < 0.4 annulus; open triangles — the galaxy 
counts in the 0.4 < (V/Y200) < 0.6 radial bin; and the solid triangles — the 0.6 < (r/r20o) < 1-0 annulus. 
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12 14 16 18 20 22 24 

R c Magnitude 

Fig. 5. — Color-magnitude diagram for A260, depicting the region where red sequence galaxies were selected. For galaxies brighter 
than R c = 18.8 (dashed vertical line), red sequence galaxies are designated as those within ±0.22 mag (±3<r) of the cluster red sequence 
relation (nearly horizontal solid line). Galaxies fainter than R c = 18.8 are selected if they lie within the region 2.5ctb-.r c rcdward of the 
red sequence and brighter than the magnitude completeness limit (solid short vertical line). 
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Fig. 6. — Composite red sequence i? c -band LF for four cluster-centric annuli: A) (r/r2oo) < 0.2, B) 0.2 < (r/r2oo) 0.4, C) 
0.4 < (r/r2oo) < 0.6, and D) 0.6 < (r/r20o) < 1-0. The short dashed line represents a Schechtcr function fit to bright-end with a fixed 
faint-end slope, a = — 1. The solid circles depict the combined net galaxy counts for all contributing clusters in each annulus. The long 
dashed line is a Schechter function fit to the faint-end, while the solid line is the sum of the two Schechter functions. 
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Fig. 7. — Superposition of the red sequence composite LF measured for four cluster-centric radial bins. The outer three LFs have been 
scaled to match the inner-most LF in the —22 < Mn a < —21 magnitude range. Plot symbols are equivalent to those defined in Figure [4] 
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Fig. 8. — Color-magnitude diagram for A260, depicting the region used to select the blue cluster galaxy population. For galaxies brighter 
than R c = 18.8 (dashed vertical line), blue galaxies are designed as those with B — R c color bluer than the lower envelope used to select 
the red sequence galaxies. Galaxies fainter than R c = 18.8 and brighter than the completeness limit (R c = 20.7, short vertical line), are 
selected if they are located in the region bluer than the lower envelope (blueward of the CMR) defining the boundary of the region used to 
select red sequence galaxies. 
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Fig. 9. — Composite blue galaxy i? c -band LF for four cluster-centric annuli: A) (r/r2oo) < 0.2, B) 0.2 < (r/r2oo) < 0.4, C) 0.4 < 
(f/1'200) < 0.6, and D) 0.6 < (r/r20o) < 1-0. The short dashed line represents a Schechter function fit to bright-end with a fixed faint-end 
slope, a = —1. The solid circles depict the combined net galaxy counts for all contributing clusters in each annulus. The long dashed line 
is a Schechter function fit to the faint-end, while the solid line is the sum of the two Schechter functions. 
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Fig. 10. — Superimposed composite blue galaxy LFs measured for four cluster-centric radial bins. The outer three LFs have been scaled 
to match the inner-most LF in the —22 < Mr c < —21 magnitude range. Plot symbols are equivalent to those defined in Figure [4] 
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Fig. 11. — Composite LF for the total, red sequence, and blue galaxy populations for four radial bins depicted in previous figures (e.g., 
see Figure [9}. The red sequence and blue LFs (represented by open squares and solid triangles, respectively) have been scaled to have the 
same net galaxy counts as the total LF (solid circles) in the —24.0 < Mr c < —16.5 magnitude range. 
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Fig. 12. — Histogram distribution of measured for 55 clusters photometrically complete to Mr c = —20 and covering a cluster-centric 
radius of (r/r20o) = 0.4. The distribution of is approximately Gaussian with (MJj } = —22.24 ± 0.06 and a dispersion of a = 0.31 

mag (dashed line). 
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Fig. 13. — Distribution of with cluster richness B gc for the 55 clusters depicted in Figure Il2l No significant correlation between 

MJj and B gc was found for our sample. 
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Fig. 14. — Mi as a function of BM-type for 54 of the 55 clusters depicted in Figure IT3l The solid line represents a least-squares fit 
to the data and indicates that MJj brightens for later BM-type. A Spearman rank-order correlation coefficient test yields that the two 
measurements are correlated at the 96% significance level. 
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Fig. 15. — Variation in M% with cluster-centric radius for the total (open squares), red sequence (solid circles), and blue galaxy 
populations (solid triangles). MJj is measured for a composite sample of 57 clusters complete to Mr c = —20. The red sequence and blue 
galaxy samples exhibit a trend in which becomes brighter with increasing cluster-centric radius. 
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Absolute R c Magnitude 

Fig. 16. — Comparison of cluster LFs with published sources: P06 - Popcsso et al. 2006, Outer — this paper, composite total LF for 
0.6 < (r/r2oo) < 10, Inner - this paper, composite total LF for (r/r2oo) < 0.2, H05 - Hansen et al. 2005, P01 - Piranomontc ct al. 2001, 
C89 - Colless 1989, D03 - De Propris ct al. 2003, G02 - Goto ct al. 2002, and B03 - Blanton ct al. 2003 (SDSS field LF; dashed line). 
The LFs have been scaled by 0.3 dex relative to each other for comparison purposes. 
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Fig. 17. — Composite total i? c -band LF for two groups of clusters containing < 10 redshift-confirmed galaxies (open squares) and > 25 
redshift-confirmed galaxies (solid circles). The LFs have been scaled to match the total composite LF in the —24 < Mr c < —17 magnitude 
range. The four radial bins are equivalent to those used in Figure [3] 
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TO OBSERVER 



Fig. 18. — Schematic diagram illustrating the geometry used to convert the projected LF to the deprojected LF. The projected central 
LF from region A will be contaminated by projected galaxies from regions Bl and B2 that lie in the cluster outskirts. The LF in region C 
can be utilized to deproject the central LF and thus minimize the influence of the contaminating galaxies. 
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Fig. 19. — Top Panel: Dcprojected total composite LFs are compared for the inner-most radial bin (r/r20o) < 0.2 (open triangles) 
with the 0.2 < (V/Y200) < 0.4 radial bin (solid triangles). The outer deprojected LF has been scaled to match the inner LF in the 
—22 < Mr c < —21 magnitude range. Middle Panel: The deprojected total composite LF (open squares) is compared to the projected 
LF (solid circles) for the (r/r20o) < 0.2 annulus. The deprojected LF has been scaled to match the projected LF in the —22 < Mr c < —21 
magnitude interval. Bottom Panel: The deprojected total composite LF (open squares) is compared to the projected LF (solid circles) 
for the 0.2 < (r/r2oo) < 0.4 annulus. The deprojected LF has been scaled to match the projected LF in the —22 < Mjj c < —21 magnitude 
range. 
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Fig. 20.— 



Same as Figure H"9l but for the red sequence composite LF. 
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Fig. 21. — Same as Figure H"9l but for the blue composite LF. 
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Fig. 22. — Vector-style color-magnitude diagram for A260. The vectors define the relative change in B — R c color and i? c -mag for the 
position of the observed galaxies to the location of the simulated galaxies. The cluster CMR is depicted by the horizontal solid line and 
the ±3c limit by the dashed line. For clarity, only a fraction of the cluster galaxies brighter than the magnitude completeness limit for this 
cluster arc displayed. 
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Fig. 23. — Composite total observed LF (solid circles) is compared to the simulated LF (open squares) for the four radial bins used 
previously: A) (r/r 200 ) < 0.2, B) 0.2 < (r/r 200 ) < 0.4, C) 0.4 < (r/r 20 o) < 0.6, and D) 0.6 < (r/r 200 ) < 1-0. The simulated LF is 
constructed by randomly changing r 2 oo by ±15%. The simulated LFs have been offset by 0.3 mag in order to assist the comparison. 
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Fig. 24. — Comparison of the total composite LF for 57 clusters that are photometrically complete to Mr c = —20 and cover the 
(r/r20o) < 0.2 annulus. The solid points represent the LF comprised by excluding the BCGs from the net galaxy counts. The inclusion of 
the BCGs in the LF is depicted by the open squares. 



